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ABSTRACT 



Over the years, there has been significant variation in the filmwise steam 
condensation data at NFS on horizontal low-integral finned tubes. With a view to 
increasing the accuracy of the data, inserts were used inside the tubes to reduce 
inside thermal resistance; however, significant discrepancies then occurred in the 
calculated outside coefficient when compared to data taken without an insert. These 
discrepancies arose due to the data reduction technique which assumes a known 
inside heat-transfer resistance and subtracts this from a measured overall resistance. 
If the assumed value on the inside is inaccurate, then the outside value is equally 
inaccurate. 

The present work uses an instrumented smooth tube to obtain accurate inside 
heat-transfer correlations both with and without inserts and uses these to obtain 
accurate outside coefficients for a family of uninstrumented finned tubes with a view 
to finding an optimum fin spacing for steam condensation. 
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I. INTRODUCTION 



A. BACKGROUND 

A reduction in size and weight of aU types of heat exchangers aboard Naval 
vessels will allow more efficient use of space. The benefits might include greater 
equipment accessibiUty for maintenance or greater heat exchanger capacity (without 
a corresponding increase in size and weight) with a corresponding increase in fuel 
efficiency. 

For the past ten years, the Naval Postgraduate School in collaboration with the 
David Taylor Research Center and the National Science Foundation has conducted 
research that is directed at the development of smaller, more efficient steam 
condensers. Improved designs can result in significant space savings, always a 
primary concern on Naval vessels, especially submarines. 

Uncertainties in past data using steam were apparently due to the lack of 
detailed information about the inside heat-transfer correlations used during the data 
reduction process. Previously, the standard Sieder-Tate correlation was assumed to 
be vah'd for the inside heat-transfer coefficient, but it may not be the best correlation 
to use with the particular test arrangement used in this research program. 

A large amount of enhanced condensation data has been collected in previous 
studies at the Naval Postgraduate School on more than 90 different condenser tubes 
of varying fin height, fin spacing, and tube material, with a view to finding an 
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optimum fin geometry for both steam and refrigerant condensation. Satisfactory 
results have been obtained with refrigerant data using R-113. However, some 
troublesome questions of possible contamination and instrument inaccuracy still 
remain with the steam data. It is felt that in order to address these questions, a 
fundamental evaluation of the heat transfer apparatus on which this data was 
collected and of the data reduction process was considered appropriate. Current 
data, taken on a carefully cleaned and caUbrated apparatus, could be compared to 
previously recorded data and a determination as to its vahdity and reproducibihty 
could be made. In addition, a thorough evaluation of the best inside heat-transfer 
correlation would lead to more reliable steam condensation results. 

B. CONDENSATION 

Condensation occurs when a vapor is cooled below its saturation temperature, 
or when a vapor/gas mixture is cooled below its dew point. Surface condensation 
occurs in condensers when a cooled surface (kept at a temperature below the 
saturation temperature of the vapor) contacts the vapor. The vapor molecules that 
contact such a surface stick to that surface and condense into liquid molecules. 
Condensation may occur in one of the following modes: filmwise, dropwise, or mixed 
mode (a combination of filmwise and dropwise) condensation. In the filmwise mode, 
the liquid wets the cold surface to form a continuous film. If the liquid does not wet 
the surface but instead forms discrete drops on the cold surface, dropwise or mixed 
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mode condensation wiU occur and is often caused by some form of contamination. 
[Ref. 1] 

The condensate forming on the tube surface offers a resistance to heat transfer 
between the vapor and the surface, which increases with the thickness of the liquid 
layer. Even though dropwise condensation results in much larger heat-transfer 
coefficients than filmwise condensation, it is difficult to maintain a stable dropwise 
condition over prolonged periods. Therefore, in most cases condenser design 
calculations are based on the assumption of filmwise condensation, resulting in lower 
heat-transfer coefficients and more conservative designs. [Ref. 1] 

In a condenser, the coolant side, tube wall, and vapor side thermal resistances 
control the heat transfer rate from vapor to coolant. Also, for experimental work 
we always use clean tubes, but in real condensers tubeside fouling can play an 
important role an increasing the coolant side resistance. The magnitudes of these 
resistances depend on the fluid, tube geometry, and flow conditions on the vapor and 
coolant side. For steam condensation, it is the coolant side thermal resistance which 
tends to dominate. Methods to lower this inside resistance include the use of inserts 
or roped tubing to promote turbulence, thereby raising the convection heat-transfer 
coefficient. However, such modifications lead to increased pressure drop through 
the tubes, which must be compensated for by providing extra pumping capacity. 
Heat transfer through the tube wall is conductive and is fixed once tube thickness 
and material are selected. The vapor side resistance is due to the condensate film 
which forms on the outside of the tube. For filmwise condensation, the outside 
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resistance can be reduced by the addition of low integral fins. These have the effect 
of not only increasing the outside surface area of the tube, but also of thinning the 
condensate film around the fins due to surface tension forces. Too small a fin 
spacing may result in condensate flooding, whereas too large a fin spacing 
approaches the smooth tube case; there should be an optimum fin spacing 
somewhere in between these two extremes. Horizontal fin spacing is therefore of 
prime importance, and finding the optimum spacing is one of the objectives of this 
long-term research program. 

C. CONDENSATION RESEARCH AT NAVAL POSTGRADUATE SCHOOL 

The research effort at NFS has included the study of differing fin dimensions 
(i.e. fin height, fin width, fin spacing) on low-integral finned horizontal tubes. 
Experimentation has included the use of three different test fluids (steam, R-113, 
and ethylene glycol) under various operational conditions using a number of different 
tube diameters. 

Van Petten [Ref. 2] provides a summary of research efforts at NPS through the 
end of 1988. Van Petten and subsequent researchers have analyzed small, medium, 
and large diameter finned tubes to find the optimum fin spacing for maximum heat- 
transfer enhancement of the fluids mentioned above. However, discrepancies found 
by Guttendorf [Ref. 3] in the data processing technique (modified Wilson plot), 
which resulted in different values of heat-transfer enhancement (for the same tube 



4 



under the same operating conditions depending on whether an insert was or was not 
used), have raised doubts about the accuracy of the inside heat-transfer correlation. 

Rouk [Ref. 4] investigated the use of an optimuation technique to predict the 
inside heat-transfer correlation. When the optimization effort proved unsuccessful, 
he next used the instrumented smooth tube data of Georgiadis [Ref. 5] to develop 
an inside heat transfer correlation, but could not find a correlation with sufficient 
accuracy based on previous data. He recommended that once an overhaul on the 
test apparatus was complete, an increase in data precision would allow the 
development of an accurate inside heat-transfer correlation. This work is a follow 
on effort to develop inside heat-transfer correlations which can predict the value of 
the inside heat-transfer coefficient with good accuracy under a variety of flow 
conditions. Once an inside correlation can be found, the object of this effort is to 
reprocess previous data and see if the discrepancies reported by former researchers 
for finned tubes on this apparatus can be rectified. 

D. OBJECTIVES 

The main objectives of this thesis were to: 

1. Disassemble and meticulously clean the apparatus to eliminate any existing 
contamination with a view to eliminating dropwise condensation problems 
experienced in the past. 

2. Carefully reassemble the apparatus using new gasket material, and make 
modifications to improve system performance. 

3. Recalibrate all system instrumentation to ensure the greatest achievable 
accuracy. 
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4. Investigate the possibility of manufacturing large, medium, and small diameter 
instrumented smooth tubes. 

5. Use the new instrumented tubes and the one existing medium diameter 
smooth instrumented tube (fabricated by Poole [Ref. 6]) to obtain accurate 
inside heat-transfer correlations for a number of insert types as well as the 
no insert condition. 

6. Evaluate the accuracy of the currently used data processing technique 
(modified Wilson plot) using instrumented tube data. 

7. Reprocess previous data using the new correlations with a view to comparing 
current and past smooth tube and finned tube data to provide continuity with 
previous studies. 
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II. LITERATURE SURVEY 



A. INTRODUCTION 

When a vapor condenses in the filmwise mode on a smooth horizontal tube it 
forms a thin continuous film of condensate on the surface of the tube. The 
condensate film thickens around the tube due to gravity. This condensate film 
provides a resistance to heat transfer which may be lowered through the use of fins. 
For quite some time, it was thought to be impractical to use finned tubes with high 
surface tension fluids such as water, due to condensate retention and flooding 
between the fins. However, a number of studies conducted on finned tubes using 
steam have shown that substantial heat-transfer enhancement may be achieved. 

A significant amount of research at the Naval Postgraduate School and 
elsewhere has addressed the issue of optimum fin height, thickness, and spacing 
required for maximum heat transfer. Yau et al [Ref. 7] reported that "with an 
increase in fin density, up to a limit (this limit is not yet known in a generahzed 
manner), the heat-transfer coefficient increases at a rate faster than the increase in 
the outside area due to the presence of fins. This additional enhancement is due to 
the thinning effect of the surface-tension forces on the condensate film. 
Unfortunately, surface-tension forces also adversely affect heat transfer by causing 
condensate to be retained between fins" [Ref. 8]. Katz et al [Ref. 9] also found that 
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on finned tubes the portion of the surface occupied by condensate is dependent upon 
the ratio of condensate surface-tension to density and the fin geometry. 

Condensate retention and the behavior of the condensate film on the tube 
surface under various conditions are critical parameters in the heat transfer process 
on horizontal finned tubes. Several models have been developed to predict this 
behavior and the reader is referred to an extensive review of horizontal finned tube 
heat transfer by Marto (Ref. 10] for a more detailed coverage of the topic. 

B. VAPOR SIDE CONSIDERATIONS 

The filmwise condensation of vapor on a horizontal tube is a complex two- 
phase heat transfer process, for which a suitably complex model would be required 
to accurately predict heat transfer performance under all conditions. 

In 1916, Nusselt [Ref. 11] set forth his theoretical work on the study of laminar 
filmwise condensation of a "stationary" vapor on a vertical or inclined plate and a 
horizontal tube. Nusselt’s simplifying assumptions included the following [Ref. 12]: 

1. Pure saturated vapor 

2. Negligible vapor velocity (U.=0) 

3. Heat transfer across the condensate film by conduction only 

4. Laminar condensate flow governed only by gravitational and viscous forces 

5. Condensate properties constant 

6. Isothermal condensing surface 

7. Negligible interface temperature drop 
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Nusselt’s result for the mean heat-transfer coefficient for a horizontal tube was 



obtained: 



=0.728 




( 2 . 1 ) 



or 



Mr =0.655 




( 2 . 2 ) 






where: 



Nu = mean Nusselt number 

kf = thermal conductivity of condensate film (W/m • k) 

Pi = condensate film density (kg/m^) 

p^ = vapor density (kg/m^) 

g = gravitational constant (9.81 m/s^) 

hfg = specific enthalpy of vaporization (J/kg) 

Pi = dynamic viscosity of condensate film (N s/m^) 

Dg = outside tube diameter (m) 

AT( = average temperature difference across condensate film (K) 
q = heat flux based on outside area (Q/A^) (W/m^) 



Many workers have improved on Nusselt’s theoretical analysis by accounting 
for some of the terms he neglected through his simplifying assumptions. However, 
equations (2.1) and (2.2) have been found to be remarkably accurate over a wide 
range of conditions for a stationary vapor. High vapor velocity can increase film 
condensation heat transfer substantially. This enhancement, which refers to the 
amount of heat transfer above or below the value predicted by the Nusselt analysis, 
is due to the effect of thinning the condensate film. However, vapor shear is the one 
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assumption which if applied can lead to significant increases in the heat-transfer 
coefficient. 

The theoretical result of Shekriladze and Gomelauri (1966) [Ref. 13], who 
considered interfacial shear stress due to vapor velocity, is shown in equation (2.3). 



=0.64(1 +(1 + 1.69f)i/2)W ^2.3) 



where: 

Nu = Nusselt number for the vapor side 
Re 24 , = two phase Reynolds number, (pj U„ DJfij) 



For steam condensation, the empirically derived correlation of Fujii et al [Ref. 
14] is shown in equation (2.4). 



Nu 



= 0.96F^5 



(2.4) 



The Nusselt expression (equation (2.1) can be expressed in similar form: 

(2.5) 

Whereas the vapor velocity, U„, cancels out in the Nusselt expression 
(stationary vapor assumption), the Fujii correlation includes the vapor velocity effect. 
Therefore we can expect equation (2.4) to more accurately predict steam-side heat 
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transfer coefficients for those cases where vapor velocity begins to have a significant 
impact. 



For further review of basic theoretical studies on the subject of laminar film 
condensation on smooth tubes the reader is referred to Rose [Ref. 15]. 

C. COOLANT SIDE CONSIDERATIONS 

For a turbulent flow regime inside a pipe (Re > 10,000), a number of coolant- 
side correlations have been used; many of these have taken the form; 



where: 

Nu = mean coolant Nusselt number for turbulent flow 
C, = correlating coefficient 
Re = coolant Reynolds number 
Pr = coolant Prandtl number 



The most common correlations with the same form as equation (2.6) are that 
of Dittus and Boelter (1930) [Ref. 16]: 



Nu=C,Rc'"Pr" 



( 2 . 6 ) 



Nu =0.023 



(2.7) 



and Colburn (1933) [Ref. 17]: 



Nu =0.023 Re^^Pr^f^ 



( 2 . 8 ) 
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Sieder and Tate (1936) [Ref. 18] applied a correction factor to equation (2.8) 
to account for cases in which the bulk to inner wall temperature difference is large 
enough to cause substantial variations in coolant viscosity as follows: 

Nu =0.021 (2.10) 

where: 

/i-c = coolant viscosity evaluated at mean bulk temperature (N-s/m^) 

IJi.^ = coolant viscosity evaluated at mean inner tube wall temperature 
(N-s/m^) 

The fluid properties in equations (2.7), (2.8), and (2.9) are evaluated at mean 
coolant bulk temperature T„^*. 

•( 7 ', » r ,)/2 ( 2 . 11 ) 

where: 

Tj = tube coolant inlet temperature (K) 

T, = tube coolant outlet temperature (K) 

Equations (2.7), (2.8) and (2.9) are vahd for Re> 10'' and 0.7 < Pr < 100, and 
were developed for long smooth pipes with no inserts [Ref. 12]. 

The use of inserts and the effect of bends close to the tube entrance region 
can affect the values of both the leading coefficient, C, , and Reynolds number 
exponent, m, in equation (2.6). One of the major focuses of this study is to 
determine the values of Q and m for the no insert, wire wrap insert, and Heatex 
insert cases. 
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Other well-known turbulent pipe flow heat-transfer correlations (i.e. Petukhov- 
Popov [Ref. 19], Sleicher-Rouse [Ref. 20], etc.) and the results of an ANL (Argonne 
National Laboratory) [Ref. 21] study which evaluated several such correlations for 
accuracy are reviewed in section VI C. 
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III. APPARATUS AND SYSTEM INSTRUMENTATION 



A. SYSTEM OVERVIEW 

The apparatus used for this research was basically the same as was used by Van 
Petten [Ref. 2] and Guttendorf [Ref. 3] with certain modifications. A system 
schematic is provided in Figure 1. Steam generated from distilled water in the .30 
m diameter pyrex glass boiler using ten 4 kW, 440 V Watlow immersion heaters was 
the working medium for this set of experiments. From the boiler section the steam 
passed up through a reducing section and a 2.13 m straight length of pyrex glass 
piping, (ID of 0.15 m), it was then turned through 180 degrees using two 90 degree 
pyrex glass elbows, and then descended down a 1.52 m straight length of pyrex glass 
piping. The steam then entered the stainless steel test section containing the 
horizontally mounted condenser tube (see Figures 1 and 2); any steam not 
condensing there was condensed in the auxiliary condenser located just beneath the 
test section. The auxiliary condenser was constructed of a single copper coil mounted 
to a stainless steel baseplate enclosed within a pyrex glass condenser section. 
Coolant flow through the auxiliary condenser was used to control system pressure, 
and all condensed liquid was returned via the condenser baseplate drain to the boiler 
section by gravity. 

Coolant for the auxiliary condenser was provided via a throttled water 
connection with associated flowmeter. Coolant flow through the single horizontal 
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Figure 1. Schematic of Single Tube Test Apparatus 
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Figure 2. Schematic of Test Section Insert 
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tube was provided by a separate system consisting of a sump tank with two 
centrifugal pumps connected in series. Coolant flow rate was measured by a 
carefully calibrated flowmeter. By varying coolant flow-rate through the single 
horizontal tube, the rate of steam condensation on the tube (and hence heat-transfer 
coefficient) could be varied. 

Non-condensible gases were removed using the vacuum pump system shown in 
Figure 3. The condensing coil for this purge system, located in the sump tank, 
served to condense steam carried through the vacuum line during the purging 
process. The vacuum line took its suction from the base of the auxiliary condenser, 
the coolest spot in the apparatus and the place where non-condensible gases (i.e. air) 
were most likely to accumulate. 

B. SYSTEM INSTRUMENTATION 

The power to the 440 V heaters was controlled through a panel mounted 
potentiometer. A description of the power calculation for input into the data 
acquisition system can be found in Poole [Ref. 6]. 

System pressure was monitored in three ways: 

1. A Setra model 204 pressure transducer 

2. A Heise solid front pressure gauge (visual reading only) 

3. System saturation pressure from vapor temperature measurement 
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Figure 3. Schematic of Purging System and Cooling Water Sump 
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DRAIN 



System vapor temperature was monitored using both teflon, and metal sheathed 
type-T copper/constantan thermocouples located juxtaposed in the test section; this 
position was just upstream of the test condenser tube. Condensate temperature was 
also monitored using a teflon coated type-T copper/constantan thermocouple located 
on the condensate return line between the auxiliary condenser and boiler. Coolant 
temperature rise in the condenser tube was measured using four methods: 

1. Two teflon coated type-T copper/constantan thermocouples 

2. Two metal sheathed type-T copper/constantan thermocouples 

3. Two Hewlett-Packard 2804A quartz crystal thermometers 

4. A ten-junction teflon coated type-T copper/constantan thermopile 

These were all placed at the inlet to and exit from the condenser tube; at the 
outlet, all thermocouples were placed just downstream of a coolant mixing chamber. 

Two data reduction programs were used to collect and reduce data on this 
apparatus; "DRPINST', and "DRPKS". The instrumented tube constructed by Poole 
[Ref. 6] was used to determine an accurate inside heat-transfer correlation for inserts 
used; this instrumented tube contained six wall thermocouples. For the instrumented 
tube the appropriate cahbration equations were accessed in the data acquisition 
program "DRPINST'. For non-instrumented tubes, the data reduction program 
"DRPKS" was used. Fluid property equations used in the data reduction programs 
are given in Appendix A and caUbrations were conducted for aU system 
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instrumentation (flowmeter, thermocouples, pressure transducer, etc.) and are 
included in Appendix B. 

The data as monitored by the aforementioned system instrumentation was 
processed by an HP-3497A data acquisition system controlled by an HP-9826A 
computer provided with the correct data acquisition program. The raw data was 
processed and stored on computer disks. Program channel assignments are given in 
Table 1. 

C. SYSTEM MODIFICATIONS 

At the beginning of this investigation the apparatus was entirely disassembled 
to facilitate complete overhaul of the system. Modifications and details of assembly 
were as follows: 



1. The apparatus was taken apart piece by piece, inspected, meticulously cleaned 
with a warm solution of Sparkleen biodegradable soap and subjected to a 
complete acetone rinse prior to reassembly. 

2. A new pyrex glass riser section above the boiler was 0.31 m shorter than the 
previous section and allowed the addition of a new aluminum stand on which 
to place the heater baseplate. This new stand allowed much easier access to 
the 440 V heater wiring plus the adjustable legs allowed level adjustment of 
the entire apparatus, to ensure proper ahgnment of the single horizontal tube. 

3. The two pyrex glass elbows were replaced. 

4. Every gasket in the system was replaced (using Buna-N rubber) and, using a 
standard star torque pattern, all flanged joints were tightened to a final 
torqued of 30 inch-pounds (manufacturer recommended maximum torque 
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Table 1. DATA ACQUISITION SYSTEM CHANNEL ASSIGNMENT 



CHANNEL 


QUANTITY 


APPLICABLE TYPE 




MEASURED 


DESIGNATION 


40 


Vapor Temperature 


T-57 

(See Appendix B.l) 


41 


Vapor Temperature 


T-56 

(See Appendix B.l) 


42 


Room Temperature 


T-58 

(See Appendix B.l) 


43 


Tube Coolant in 1 


T-58 


44 


Tube Coolant out 1 


T-58 


45 


Tube Coolant in 2 


T-55 

(See Appendix B.l) 


46 


Tube Coolant out 2 


T-55 


47 


Condensate Return 


T-58 


48 


Instrumented Tube 


T-57 or T-55 
as applicable 


49 


Instrumented Tube 


T-57 or T-55 
as applicable 


50 


Instrumented Tube 


T-57 or T-55 
as applicable 


51 


Instrumented Tube 


T-57 or T-55 
as applicable 


52 


Instrumented Tube 


T-57 or T-55 
as applicable 


53 


Instrumented Tube 


T-57 or T-55 
as applicable 


60 


10 - Junction 


T-61 




Thermopile 


(See Appendbc B.l) 


61 


Voltage 


NA 


62 


Current 


NA 


64 


Pressure Transducer 


(See Appendix B.4) 
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specification was 60 inch-pounds). The bolts holding the flanges together 
should be checked for tightness on a regular basis, since the thermal cycling 
of the apparatus has been shown to cause loosening of nut and bolt 
assemblies. 

5. The previous vacuum pump was replaced with a Cast model 2567-V108 
vacuum pump which could draw vacuum to 130 mmHg with an installed check 
valve to prevent pump back spin when the pump was stopped. The new 
pump could draw vacuum much more rapidly, but could not operate at less 
than 130 mmHg. Remaining non-condensible gases were removed by flushing 
the system with steam. The steam flushing procedure for removal of non- 
condensible gases is given in the operating procedures section. Once the 
system is completely filled with steam, operating pressures well below 130 
mmHg could be achieved utilizing the auxiliary condenser. 

6. The double coil auxiliary condenser was replaced with a single coil used 
originally by Van Petten [Ref. 2]. The single coil was not coated with the 
special oxide coating used by Guttendorf [Ref. 3]. It was felt to be superior 
to the double coil in that the baseplate welds were of much higher quality and 
were preferred on the basis of vacuum tightness. 

7. The aluminum side plates attached to the pyrex glass auxiliary condenser 
housing were replaced with new stainless steel side plates with penetrations 
for pressure bleed, vacuum line, and a pressure transducer. These three 
penetrations were fitted with screw threaded stainless steel connectors. The 
stainless steel connectors were heli-arc welded in place. Prior to this 
modification (completed 24 January 1992) a leak test conducted from 21 
December 1991 to 2 January 1992 revealed a mean vacuum leak rate of ~ 3.4 
mmHg per day (see Figure C.l). A subsequent leak test conducted from 6 
February 1992 to 19 February 1992 showed an improvement in the mean leak 
rate to —1.7 mmHg per day (see Figure C.2). 

8. System instrument modifications included the addition of the Setra pressure 
transducer and the Heise pressure gauge, and the removal of the mercury 
manometer. All system instrumentation was recah'brated and the results 
incorporated into the data reduction programs. 

9. Finally the apparatus was lagged with Halstead insulating foam to reduce heat 
loss as much as possible. The test section, which was left uncovered 
previously, was also lagged. 
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IV EQUIPMENT OPERATION AND EXPERIMENTAL PROCEDURE 



A. SYSTEM STARTUP AND SHUTDOWN PROCEDURES 

Startup of the system is accomplished in the following manner: 

1. Ensure distilled water level in the boiler is 4 to 6 inches above the top of the 
heating elements. The boiler is filled by gravity drain via a hose connection 
from the distilled water tank to the boiler fill valve. Ensure the vent valve on 
the side of the auxiliary condenser is open when filling or draining the boiler. 
The boiler may be drained by removing the hose connection and opening the 
fill valve, which allows drainage into the trench directly beneath the boiler. 

2. Once the boiler is filled to the appropriate level, shut the boiler fill valve and 
the distilled water tank valve. 

3. Shut the system vent valve. 

4. Turn on the data acquisition system, computer and printer. Load the 
appropriate program (DRPINST or DRPKS) and check for proper operation. 
TTien check all thermocouple outputs, by stepping through the appropriate 
data acquisition system channels, to verify that all are registering ambient 
temperature. 

5. Open the fill valve to the coolant water sump tank to a level such that the 
tank overflow drain box does not overflow (the valve is located between the 
boiler control panel and heat pipe apparatus). 

6. Turn on the coohng water supply pumps and adjust the tube flow rate from 
20% to 60% of the rotameter setting and check for leaks. Reset flow rate to 
desired level. 

7. Open valves from tap water system to auxiliary condenser and adjust coolant 
flow rate to at least 30% and check for leaks. Reset flow rate to at least 
10 %. 

8. Energize heaters and adjust voltage to approximately 50 volts (40 volts if the 
system is already at vacuum below 100 mmHg to limit the vibrational shock 
to the system from oversized vapor bubble formation). To energize the 
heaters there are three switches which must be placed in the on position. The 
first is located in power panel p5 located in the main hallway adjacent to the 
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lab and is labeled switch 3 / heater controller room 106. The second is the 
heater load bank circuit breaker located on the side of the boiler control 
panel. The third is the condensing rig boiler power switch on the front of the 
boiler control panel. Increase the voltage gradually in 10 volt increments to 
the desired level. 

9. Turn on the vacuum pump and open the vacuum line valve. Allow the 
vacuum pump to run until system pressure is below 3 psi, then shut the 
vacuum line valve just prior to turning off the vacuum pump. 

10. As system warmup continues and pressure increases to above 4 psi energize 
the vacuum pump as necessary to flush the non-condensible gases out of the 
system through the vacuum line by forcing the gases out with steam. To 
ensure that the non-condensibles gather at the base of the auxiliary 
condenser, where the vacuum line suction is located, ensure that the 
horizontal tube is not supplied with coolant flow, and adjust coolant flow 
through the auxiliary condenser as necessary to ensure steam is filhng the 
entire system. The auxiliary condenser may be touched lightly by hand, along 
its entire length, to ensure the system is completely filled with steam; any cool 
spots indicate the presence of non-condensible gases which means that the 
flushing process is not complete. The flushing process takes 15 to 30 minutes 
to accomplish, and should be repeated periodically for long periods of 
operation. 

11. At the conclusion of the flushing process, shut the vacuum line valve and 
secure the vacuum pump. 

12. In order to ensure that filmwise condensation occurs on the tube, coolant flow 
through the tube must be initiated as follows: 



a. Allow the apparatus vapor temperature measurement (channel 40) to 
reach at least 3800 microvolts. 

b. Cut in the auxiliary coolant flow (50% or 60% level) to cool the vapor 
temperature to roughly 3200 microvolts. 

c. Secure coolant flow through the auxiliary condenser, and allow the 
vapor temperature level to climb to about 3700-3800 microvolts, which 
allows a steam blanket to cover the tube. 

d. Initiate coolant flow through the single horizontal tube at the 80% level. 
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e. Cut in coolant flow to the auxiliary condenser to control pressure, and 
observe the condensation process to ensure that a condensate film has 
formed on the tube. 

13. Run software program DRPINST for an instrumented tube (DRPKS for 
uninstrumented tubes) by pressing "run" on the keyboard. 

To take data for an instrumented tube, the questions for DRPINST can be 

answered as follows: 

• Select fluid ... Enter 0 for steam 

• Select option ... Enter 1 to take new data 

• Enter month, date and time ... Press enter 

• Enter input mode ... Enter 1 for new data 

• Give a name for the raw data file ... Enter name 

• Enter geometry code ... Enter 1 for finned, 0 for plain 

• Select insert type ... Enter 0 for none, 1 for twisted tape, 2 for wire wrap, 3 for 

Heatex 

• No. of thermocouples in wall? ... Enter 4, 5, or 6 depending on the tube 

• Select tube diameter type ... Enter 2 for medium 

• Enter pressure condition ... Enter 0 for vacuum, 1 for atmospheric 

• Give a name for the wall temperature file ... Enter name 

• Select input ... Enter 1 for short, 2 for long, or 3 for raw data 

• Like to check NG (non-condensible gas) concentration ... Enter 1 for yes, 2 for 

no; you must answer yes for the first data point 

• Enter flowmeter reading (%) ... Enter 2 digit number (i.e. 20 or 58 etc.) 

• Connect voltage line ... Flip the voltage line toggle switch, located on the power 

control panel, to the on position and press enter 
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• Disconnect voltage line ... Flip the voltage line toggle switch off and press enter 

• Enter pressure gauge reading (Pga) ... Enter reading off gauge in psi 

• Select measurement ... Enter 0 for teflon, 1 for metal sheath, 2 for quartz, 

3 for thermopile 

• Change TCOOL rise? ... Enter 1 for yes, 2 for no 

• OK to store this data set? ... Enter 1 for yes, 0 for no 

• Will there be another run? ... Enter 1 for yes, 0 for no; starts at check NG 

concentration for following runs. 

To take data for an uninstrumented tube the questions for DRPKS can be 
answered as follows: 

• Select fluid ... Enter 0 for steam 

• Select option ... Enter 0 to take new data 

• Enter month, date and time ... Press enter 

• Enter disk number ... Enter number 

• Enter input mode ... Enter 0 for new data 

• Select Ci ... Enter 0 to find a Ci value, 2 to use a Ci value stored in the 

program. 

• Give a name for the raw data file ... Enter name 

• Enter geometry code ... Enter 1 for finned, 0 for plain 

• Enter insert type ... Enter 0 for none, 1 for twisted type, 2 for wire wrap, 3 for 

Heatex 

• Select tube type ... Enter 0 for thick wall (only thick wall tubes were tested) 

• Select material code ... Enter 0 for copper (only copper tubes were tested) 
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• Select tube diameter type ... Enter 1 for medium (no small or large diameter 

tubes were tested) 

• Enter pressure condition ... Enter 0 for vacuum, 1 for atmospheric 

• Want to create a file for NR vs F? ... Enter 1 for yes, 0 for no 

• Give a name for plot data file ... Enter name; easiest to use the raw data file 

name preceded by a P 

• Select output ... Enter 0 for short, 1 for long, 2 for raw data 

• Like to check NG concentration ... Enter 1 for yes, 2 for no; you must answer 

yes for the first data point 

• Enter flowmeter reading (%) ... Enter 2 digit number (i.e. 20 or 60 etc.) 

• Connect voltage line ... Flip the voltage line toggle switch on and press enter 

• Disconnect voltage line ... Flip the voltage line toggle switch off and press enter 

• Enter pressure gauge reading (Pga) ... Enter reading off gauge in psi 

• Select measurement ... Enter 0 for teflon, 1 for metal sheath, 2 for quartz, 

3 for thermopile 

• Change TCOOL rise ... Enter 1 for yes, 2 for no 

• OK to store this data set ... Enter 1 for yes, 0 for no 

• Will there be another run ... Enter 1 for yes, 0 for no; starts at check NG 

concentration for following runs. 

14. Only answer the program questions up to "Enter flowmeter readings". 
Monitor system temperature using the vapor thermocouple voltage reading 
(the program automatically resets to channel 40) closely until system warmup 
is complete. 

15. Monitor system temperature and pressure carefully to prevent a system over 
pressure during warmup (especially at atmospheric conditions). 
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16. If conducting a vacuum run, gradually adjust voltage to 90 volts (usually in 10 
volt increments). Obtain the desired operating condition by manually 
controlling coolant flow through the auxiliary condenser until channel 40 
reads 1970 ± 20 microvolts (--48'’C). Vapor velocity ~ 2 m/s. 

17. If conducting an atmospheric run, gradually adjust voltage to 175 volts from 
the 90 volt level in 10-20 volt increments. Again the desired operating 
condition is obtained by manually controlling coolant flow through the 
auxiliary condenser until channel 40 reads 4280 ± 20 microvolts (~ 100°C). 
Vapor velocity ~ 1 m/s. 

18. Monitor the condensation process using the glass viewing window periodically 
to ensure that filmwise condensation is maintained. To clear the viewing 
window of fog and moisture increase coolant flow through the auxiliary 
condenser briefly to 50% or 60%, then reset to desired flow rate. 

19. When taking readings be sure to check the flowmeter setting prior to entering 
it into the computer (it has a tendency to fluctuate slightly). 

20. If conducting vacuum and atmospheric runs on the same day always conduct 
the vacuum run first. If the atmospheric run is done first it takes too long for 
the system to cool down to vacuum operating temperatures. 

The system is secured in the following manner: 

1. Secure power to the heating elements. 

2. Secure coolant flow through the tube, through the auxiliary condenser, and 
to the sump tank. 

3. If desired to maintain the system at vacuum conditions until the next run the 
shutdown is complete. Continued cooling water circulation may be used to 
assist in cooling down the system. 

4. To bring the system back to atmospheric conditions slowly open the vent 
valve. 

5. The data acquisition system may be turned off whenever it is not necessary 
to monitor system parameters. 

6. Periodically change distilled water in the boiler. 
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7. If an emergency should arise such as abrupt overpressurization or breakage, 
immediately secure power to the heaters and open the vent valve, then let the 
system cool down before checking the apparatus for damage. 

B. EXPERIMENTAL PROCEDURES AND OBSERVATIONS 

Water is a poor wetting medium and therefore great care was taken to ensure 

that uniform filmwise condensation was the only condensation mode occurring during 

a data run. Even though the apparatus was meticulously cleaned (as mentioned 

previously), a continuing problem with dropwise condensation manifested itself. 

Subsequent to steam cleaning the system with a Sparkleen soap solution, by 

operating the system with a soapy solution in the boiler (the solution bubbled 

through the entire apparatus), dropwise condensation was observed on the installed 

instrumented tube. After taking some data when in the dropwise condition, the tube 

was removed and rigorously cleaned using a warm Alconox soap solution with a 

scrub brush. However, after observing the filmwise mode initially, the condensation 

mechanism soon transitioned to mixed mode and then back to the dropwise mode. 

Since only a filmwise condition over several hours would suffice, the tube 

chemical treatment procedure used by Guttendorf [Ref. 3] and several other 

researchers was used to produce filmwise condensation. The tube was chemically 

treated prior to installation as follows: 

1. Clean the internal and especially the external surfaces of the tube using a soft 
brush and mild soap (using the Alconox detergent in warm water), rinse with 
acetone then rinse thoroughly with distilled water. Repeat the cleaning 
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procedure until the distilled water rinse perfectly wets the tube surface; any 
breaks in the wetting film at this point are likely to result in dropwise 
condensation spots once the tube is installed in the apparatus. 

2. Place the tube in a steam bath. 

3. Mix equal amounts of ethyl alcohol and a 50% by weight solution of sodium 
hydroxide. Keep the solution warm so that a watery consistency is 
maintained. 

4. Apply the solution to the tube with a small paint brush, retaining the tube in 
the steam bath. If the tube has not been treated previously, apply a coating 
of the solution every 10 minutes for an hour. If the tube has been previously 
treated, apply a coating every 5 minutes for a period of 20 minutes. 

5. Remove the tube from the steam bath and thoroughly rinse the tube with 
distilled water to remove any excess solution. Install the tube in the test 
section immediately, being careful not to touch the tube surface. Oil or dirt 
from any source may contaminate the tube surface and result in mixed mode 
or dropwise condensation. 



The oxide layer which forms on the tube is very thin, and has negligible thermal 
resistance and high wetting characteristics. 

Once the tube was installed in the apparatus (with the desired insert in place), 
the system startup procedure outlined in section IV A was followed to take data at 
desired conditions. 

At vacuum conditions, when single tube coolant flow was initiated with vapor 
velocity at ~ 2 m/s, the condensation on the tube did not develop as a perfect film 
but instead left patches where the film was broken. These patches, or streaks 
seemed to occur at regular intervals, and it was postulated that they were due to 
vortex shedding of vapor around the tube. Therefore, the procedure in the startup 
section IV A, step 12, was used to promote the development of uniform filmwise 
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condensation by inducing a stationary vapor condition around the tube. This allowed 
a steam blanket to form around the tube prior to coolant flow initiation. After flow 
initiation the appearance of the condensate film on the tube surface was continuous 
with no breaks. Momentary film instabilities were observed at vacuum conditions at 
pressures below ~ 20 kPa after a continuous film was established at higher pressures. 
These instabilities seemed to interrupt the film sheet only for an instant and then 
disappear, and may have been caused by vortex shedding of the vapor around the 
tube as already mentioned. A possible mechanism to explain these instabihties is 
that the higher vapor velocity at vacuum conditions momentarily thins the condensate 
film via vortex shedding, yet this thinning effect is overcome by surface tension forces 
in the film sheet which tend to restore the continuous film. These instabilities could 
only be seen for an instant and then would vanish, being very transitory in nature. 
There also did not appear to be any pattern whatsoever to the instability formation. 
The instabihties were not observed at pressures above ~ 20 kPa. 

The data taking regimen for each data set involved starting, then verifying the 
existence of a filmwise condensation condition, then taking data at flow rates (in %) 
of 80, 70, 60, 50, 40, 30, and 20 then back to 80 and 50 to check for repeatability 
within the data set. Two data points were taken at each of the first seven data points 
and one each for the last two, which gave a total of 16 data points. It was usually 
quite clear from the two comparison points, and from data taken previously under 
similar conditions, whether the data set should be rejected or accepted. After tube 
installation, the appearance of one or more small patches (breaks in the film) after 
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— 7-10 hours of operation signaled the beginning of tube contamination which got 
worse with time. The tube would then be removed and cleaned. 

C. TUBES TESTED 

The data taken during this study involved extensive use of the instrumented 
smooth tube (SOI) fabricated by Poole [Ref. 6], (six wall thermocouples spaced 60 
degrees apart placed at midwall and midlength). Due to excessive thermocouple 
wear, only 5 thermocouples in this instrumented tube functioned properly. The tube 
was positioned in the apparatus to make optimum use of functioning thermocouples. 
The preferred arrangement placed 4 thermocouples at 10°, 190°, 250°, and 310° from 
the top dead center position of the tube. This arrangement provided readings from 
the top and bottom of the tube and two intermediate points, giving the most accurate 
mean tube wall temperature and best temperature profile readings. At the 
conclusion of this study, the manufacture of the new instrumented tubes was not yet 
complete; progress to date (March 1992) is recorded in Appendix F. 

Data was also taken on a uninstrumented smooth tube (S02) and four finned 
tubes (SOS, S04, SOS, S06) with fin spacings of 0.5 mm, 1.0 mm, and 1.5 mm, and 2.0 
mm. All tubes tested were classified as medium tubes, with an outside diameter of 
19.05 mm, and an inside diameter of 12.70 mm. The four finned tubes all had the 
same fin height of 1.0 mm, and fin width of 1.0 mm. Data runs were taken either 
with no insert, with the wire wrap or the Heatex insert installed.By more efficiently 
mixing the coolant an insert significantly increases the inside heat-transfer coefficient. 
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The Heatex insert consists of a central wire core onto which are wound a series 



of wire loops, each inclined at a common angle to the core. The loops come into 
direct contact with the tube wall, and each loop provides a significant amount of 
coolant mixing as the coolant flows through the loop mesh. [Ref. 12]. 

The wire wrap insert was a copper wire spirally wrapped around a central 
stainless steel rod with a uniform pitch. This insert induced a swirling coolant 
motion, which enhanced turbulent mixing within the tube. This particular wire wrap 
insert was the same insert used by Guttendorf [Ref. 3], Coumbes [Ref. 22], and Van 
Petten [Ref. 2] to collect data on the medium family of tubes. 

A summary of data runs is given in Table 2. 



Table 2. DATA RUN SUMMARY 



Tube 



Filename 



Pressure 

(kPa) 



Vapor 

Velocity 



Insert 



SOI 



DSOIMVSH2 

DSOIMASH4 

DSOIMASH5 

FSOIMASH3 

FIMAVSHl 

FIMAVSH2 

FIMAVSH3 

FIMAVSH4 

FIMAVSH5 

FIMAVSH6 

FIMAVSH7 

FIMASW3 

FIMASW4 



11 

101 

101 

11 

28 

41 

69 

101 

28 

41 

69 

101 

101 



(m/s) 

2.17 

1.06 

1.06 

2.11 

3.53 
2.45 

1.53 
1.08 
1.01 
1.04 
1.02 
1.07 
1.07 



wire wrap 
wire wrap 



Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 
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Tube 


Filename 


Pressure 

(kPa) 


Vapor 

Velocity 

(m/s) 


Insert 




FIMASW5 


101 


1.06 


wire wrap 




COMPAIWl 


101 


1.07 


wire wrap 




COMPAIW2 


101 


1.07 


wire wrap 




FIMASN4 


101 


1.08 


none 




FIMASN5 


101 


1.07 


none 




FIMASN6 


101 


1.08 


none 




COMP AIN 1 


101 


1.08 


none 




COMPAIN2 


101 


1.08 


none 


S02 


FNMAVSHl 


15 


6.24 


Heatex 




FNMAVSH2 


41 


2.42 


Heatex 




FNMAVSH3 


69 


1.51 


Heatex 




FNMAVSH4 


102 


1.03 


Heatex 




FNMAVSH5 


28 


1.02 


Heatex 




FNMAVSH6 


41 


1.03 


Heatex 




FNMAVSH7 


69 


1.02 


Heatex 




FNMAVSH8 


101 


1.07 


Heatex 


S03 


FIMAF051 


101 


1.06 


Heatex 




FIMAF052 


101 


1.06 


Heatex 


S04 


FIMAFlOl 


101 


1.07 


Heatex 




FIMAF102 


101 


1.06 


Heatex 


SOS 


FIMAF151 


101 


1.07 


Heatex 




FIMAF152 


101 


1.07 


Heatex 


S06 


FIMAF201 


101 


1.07 


Heatex 




FIMAF202 


101 


1.07 


Heatex 
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Tube descriptions: 

501 instrumented smooth tube fabricated by Poole 

502 uninstrumented smooth tube 

SOS 0.5 mm finned tube 

504 1.0 mm finned tube 

505 1.5 mm finned tube (this tube designated as F096 by Guttendorf; also a 
similar 1.5 mm finned tube designated as F006 by Guttendorf and Van 
Petten was not tested); the 1.5 mm finned tube results were not included 
in subsequent finned tube analysis since the reason for differences in 
experimental results for these two tubes has not yet been resolved. 

506 2.0 mm finned tube 
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V. THEORETICAL BACKGROUND AND DATA REDUCTION PROCEDURES 



A. THEORETICAL BACKGROUND 

The overall or total resistance to heat transfer from vapor to coolant consists 
of the sum of the vapor side resistance (Rq), the tube wall resistance (R*), and the 
coolant side resistance (R,); this neglects any fouling resistance since clean tubes are 
always used. 



The vapor and coolant side resistances are convective in nature and may be 
expressed by the reciprocal of their respective heat-transfer coefficient and surface 
area product. 




(5.2) 




(5.3) 



where; 

outside vapor side resistance to heat transfer (K/W) 
outside heat- transfer coefficient (W/m^ K) 
effective outside surface area (m^) 
inside coolant side resistance to heat transfer (KAV) 
inside heat-transfer coefficient (W/m^ K) 
effective inside surface area (m^) 

The tube wall resistance is conductive in nature and is represented by the radial 

conduction equation. 
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(5.4) 



in^ 

D. 



2nLk 



m 



where: 

= tube wall resistance (K/W) 

D, = outside or root (if finned) diameter (m) 

D, = inside tube diameter (m) 

L = active condensing length (133 mm) 

= thermal conductivity of tube wall (W/m-K) 

The effective outside area of the tube is calculated using the following 

expression: 



A^^nD^L (5.5) 

where: 

Df = outside or root (if finned) diameter (m) 

L = active condensing length (133 mm) 

The effective inside area includes the inside surface area involving the active 
condensing length and the inside surface area of the insulated inlet and outlet 
portions of the tube. These portions of the tube act as fins and remove heat via 
axial conduction. The extended fin assumption with associated fin efficiencies was 
used to account for these end losses. 



yl. = 7tD.(I + Ti^I^+Ti2l2) 



(5.6) 



where: 

D, = inside diameter of tube (m) 

L = active tube condensing length (m) 

T/i = fin efficiency of inlet portion of tube 
Lj = length of inlet portion of tube (m) 

7)2 = fin efficiency of outlet portion of tube 

L 2 = length of outlet portion of tube (m) 

The overall thermal resistance to heat transfer may be expressed by: 
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1 



(5.7) 



o o 




Substituting equations (5.2) and (5.3) into equation (5.7) yields: 




(5.8) 



where: 



Uo = overall heat-transfer coefficient (W/m^ K) 

Ao = effective outside surface area (nr) 
ho = outside heat-transfer coefficient (W/m^-K) 

= Tube wall resistance (K/W) 
h, = inside heat-transfer coefficient (W/m^ K) 

The single tube condenser apparatus uses the log mean temperature difference 
(LMTD) analysis for calculation of the heat transfer between the hot vapor and cold 



Q = heat transfer rate to the cooling water (W) 

Uo = overall heat-transfer coefficient (W/m^ k) 

Ao = effective outside surface area (mr) 

LMTD = log mean temperature difference between vapor and coolant (K) 
The log mean temperature difference (LMTD) is given by: 



coolant. 



Q = U^A^{LMTD) 



(5.9) 



where: 



LMTD = 



In 




(5.10) 
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where: 



T, = coolant inlet temperature (K) 

Tj = coolant outlet temperature (K) 
Ts 3 , = vapor saturation temperature (K) 



In this and previous studies at NPS, the quartz thermometer output for Tj and 
T 2 were used in the calculations for the coolant temperature rise, and the saturation 
temperature, T^g,, was measured using the vapor thermocouple (channel 40). 

The total heat transfer across the tube is experimentally determined by 
measuring the mass flow rate of fluid through the tube and its accompanying 
temperature rise. 



(3=mc/rj-r,) (5.11) 

where: 

Q = heat transfer rate (W) 
m = mass flow rate of coolant (kg/s) 

Cp = specific heat of coolant at constant pressure (J/kg-K) 

T, = coolant inlet temperature (K) 

T 2 = coolant outlet temperature (K) 

Equation (5.11) may be used directly from the experimental data. The 
resultant heat transfer rate, Q, is then substituted into equation (5.9) to find the 
overall heat-transfer coefficient, U^. 




Q 

A^(LMTD) 



(5.12) 



where: 

Q = heat transfer rate from eq. (5.11) (W) 

A„ = effective outside surface area (m^) 

LM I'D = log mean temperature difference; eq. (5.10) (K) 
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Since Uq, A„, and A, are known quantities, this leaves only two unknowns 
in equation (5.8), and h„ the outside and inside heat transfer coefficients. 

Often the coolant side thermal resistance is dominant, and inserts such as those 
mentioned in section IV C are used to lower the inside resistance. This allows a 
more accurate computation of the outside heat-transfer coefficient, h^,, when using 
the modified Wilson plot technique mentioned in section V B. Vapor side heat 
transfer may also be enhanced through the use of fins, drainage strips, or dropwise 
condensation promoters. 

B. MODIFIED WILSON PLOT TECHNIQUE 

The ideal way to solve for h„ and h, in equation (5.8) is through the use of 
instrumented, which accurately determine a mean tube wall temperature. The inside 
and outside mean tube wall temperatures may then be obtained directly by assuming 
a linear temperature profile across the wall. Since the vapor temperature and mean 
coolant temperature are known, the inside and outside heat-transfer coefficients may 
then be calculated directly using equation (5.13). 

q = hLT (513) 

where: 

q = heat flux (W/m^) 

h = heat-transfer coefficient (h, for inside, h^ for outside) (W/m^ K) 

AT = temperature difference across resistive medium (AT = Tsa,-T^au^„u^ijj 
for the vapor side, and AT = T,^, 3 U j„si(,e - T^o, 3 „, for the coolant side) 

For data collection on a large number of tubes, the use of instrumented tubes 
is impractical due to the high cost and difficulty involved in manufacturing so, many 
tubes. Therefore, the modified Wilson plot technique was developed, which solves 
for the inside and outside heat-transfer coefficients simultaneously without using wall 
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thermocouples. To obtain the most accurate results with this method, it is necessary 
that the inside and outside coefficients be relatively equal in magnitude. 

The modified Wilson plot technique requires that the "form" of the equation 
for both the inside and outside heat-transfer coefficients be known. The Nusselt 
theory and Sieder-Tate correlation are used to represent the "form" of the outside 
and inside heat-transfer coefficients respectively. The Nusselt theory when based on 
q can be represented by: 



h^ = a 




(5.14) 



where: 

h„ = outside heat-transfer coefficient (W/m^ K); based on q 
a = dimensionless coefficient 

kf = thermal conductivity of condensate film (W/m-K) 

Pf = condensate film density (kg/m^) 

g = gravitational constant (9.81 m/s^) 

hfg = specific enthalpy of vaporization (J/kg) 

= dynamic viscosity of condensate film (N s/m^) 

Do = outside tube diameter (m) 

q = heat flux based on outside area (Q/Aq) (W/m^) 

The Sieder-Tate correlation may be represented by: 









\ 0.14 






= C.Q 



(5.15) 



where: 

h, = inside heat transfer coefficient (W/m^ K) 

C, = Sieder-Tate leading coefficient 

kj = thermal conductivity of coohng water (W/m-K) 

D, = inside tube diameter (m) 

Re = Reynolds number 

Pr = Prandtl number 

p-o = dynamic viscosity of coohng water at bulk temperature (N-s/m^) 
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= dynamic viscosity of cooling water at mean inner wall temperature 
(Ns/m^) 

Substituting equations (5.14) and (5.15) into equation (5.8) gives: 



U. 






z= 



1 

+ — 
CjQi4. a 



(5.16) 



By letting: 



y= 






(5.17) 



X = 



Afi 



(5.18) 




(5.19) 



b = 



1 



a 



then equation (5.16) reduces to: 



Y=mX+b 



(5.20) 



(5.21) 



The parameters H and Z are both temperature dependent, and must be 
determined iteratively. A least-squares fit of equation (5.21) is used to determine 
Q and a. Once C, is known, h, can be calculated using equation (5.15). With h, and 



42 



Uq known, the value of \ can then be easily determined by rearranging equation 
(5.7) as in equation (5.22), or by using equation (5.14) with the value of a known. 




(5,22) 



C. INSTRUMENTED TUBE IMPROVEMENTS FOR DATA REDUCTION 

In previous work at NPS, the standard form of the Sieder-Tate equation was 
used with a Reynolds number exponent of 0.8, as in equation (2.10). One of the 
aims of this thesis was to use an instrumented tube to directly determine the inside 
and outside coefficients, h; and h<> , and then use the data to determine a more 
"exact" form of the Sieder-Tate-type equation to be used for each insert. The 
coolant side correlations mentioned in section II C were based on a long, straight 
entrance length. The sharp 90° bend just prior to the test section tube entrance 
undoubtedly creates entrance effects which lead to discrepancies between our 
experimental data and heat transfer behavior predicted by the Sieder-Tate 
correlation. 

Assuming the final form of the inside heat-transfer correlation to be that of 
equation (2.10) gives the following: 



Nu=C.Re'^pM— 



f ^ \ 0.14 



(5.23) 




where: 



m 



= Reynolds number exponent to be determined 
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Rearranging equation (5.23) gives: 



Nu 



Prin 






\0.14 



= C.Re”* 






(5.24) 



Taking the natural log of equation (5.24) yields: 



In 



/ 


Nu 


\ 




u ) 


0.14 


Prift 


Pc 




\ 




/ 



= \n a + m ]n Re 



(5.25) 



Equation (5.25) is in the form of a linear equation, and by plotting In 
(Nu/(Pr^^(/xy/i.^)®’‘‘*)) versus In Re, the slope and intercept, namely the Reynolds 
number exponent and Sieder-Tate coefficient may be determined from the 
instrumented data. 

With the unknown parameters of equation (5.23) determined, the new inside 
heat-transfer correlations (one for each insert) could be used in the data reduction 
program to give the value of the inside heat-transfer coefficient, h„ directly and to 
provide a more accurate calculation of the outside heat-transfer coefficient h^. 

D. ENHANCEMENT RATIO 

Following the development of Van Petten [Ref.2], and Nusselt theory, 
experimental data can be curve fitted, using a least-squares analysis, to an equation 
of the following form: 
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q = a A r” 



(5.26) 



where: 

q = heat flux based on outside area (Q/A„) (W/m^ 
AT = temperature drop across the condensate film (IC) 

Substituting equation (5.12) into his expression yields: 



O 



(5.27) 



From Nusselt theory n = 0.75, therefore the enhancement ratio, based on 
constant temperature drop across the condensate film, can be expressed as: 



€ 



^T 




a 



1 



a. 



(5.28) 



where: 

= enhancement ratio based on constant temperature drop across the 
condensate film 

/ = subscript denoting finned tube 

s = subscript denoting smooth tube 

ho = outside heat-transfer coefficient (W/m^-K) 

a = constant of proportionality introduced in equation (5.26) 



Also for constant AT; using equation (5.14): 











1/3 










1 







(5.29) 



If the heat flux is kept constant the values of Zj and remain equal, which 
results in equation (5.30). 
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(5.30) 



€ = ^1 - “/^/ - “/ 
^ h ^ a Z a 

os s s s 



Combining equations (5.29 and (5.30) gives the relationship between and 
in equation (5.33). 



6 



f\ 



a. 



a 

s 



a 



1 



14/3 



a. 






Note that e^j and are independent of q and AT. 



(5.31) 



(5.32) 



(5.33) 
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VI. RESULTS AND DISCUSSION 



A. DROPWISE CONDENSATION 

As mentioned in section IV B, a dropwise condensation condition was obtained 
initially. Data were taken during dropwise conditions at vacuum and atmospheric 
pressure with the instrumented tube, and this data is compared to filmwise data in 
Figure 4. The figure shows a marked contrast between filmwise and dropwise 
condensation data. The dropwise heat-transfer enhancement, compared to the 
filmwise data, varied from ~2 to 7 for vacuum conditions and —9 to 10 for 
atmospheric conditions. Marto et al [Ref. 23] studied the use of organic coatings for 
the promotion of dropwise condensation of steam, and obtained an outside heat- 
transfer coefficient of —55 kW/m^-K for a Fluoroacrylic coating compared to —30 to 
85 kW/m^ K found under vacuum conditions in this study (both taken at P— 11 kPa 
and vapor velocity —2 m/s). The difficulty in accurately measuring the temperature 
drop across the condensate film for the dropwise condition is illustrated by the large 
amount of scatter in the dropwise data; this problem is caused by the lack of a stable 
film and the intermittent presence of drops near the instrumented tube wall 
thermocouples. The exact cause of the dropwise conditions was never determined. 
Clearly, some organic contamination either from the boiler feed water or from the 
gasket material was depositing on the test tube. Since this thesis was devoted to 
filmwise condensation, great efforts were made to clean the test tube and prepare 
it chemically so that the condensate would wet the surface. 
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Figure 4. Comparison of Dropwise and Filmwise Condensation Data (Smooth 
Instrumented Tube, Heatex Insert) 
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B. INSTRUMENTED TUBE RESULTS 



After resorting to the oxide coating procedure mentioned previously to mitigate 
the dropwise contamination problem, a series of runs were made with the 
instrumented tube fabricated by Poole [Ref. 6] with the Heatex insert installed. 
Using a mean wall temperature, the inside and outside heat-transfer coefficients 
could be evaluated directly. 

Figure 5 shows filmwise condensation data taken at various pressures and vapor 
velocities. The four data sets for increasing vapor velocity (from 1 to 3.5 m/s) and 
decreasing vapor pressure (from 101 to 28 kPa) show the effect of vapor shear 
thinning the condensate film, giving an increase in the outside heat-transfer 
coefficient. This series of data runs was taken by maintaining the heater voltage at 
175 volts and adjusting coolant flow through the auxiliary condenser to control 
pressure. 

The four data sets in Figure 5 with constant vapor velocity (~1 m/s) and 
decreasing vapor pressure (from 101 to 28 kPa) show reduction of the outside heat- 
transfer coefficient with decreasing saturation temperature (due to decreasing 
saturation pressure). This effect is thought to be due to an increase in condensate 
viscosity at lower temperatures which tends to prevent the condensate from flowing 
around and draining from the tube as easily as at higher temperatures. The resulting 
condensate film thickening provides an additional resistance to heat transfer thereby 
lowering the outside heat-transfer coefficient. This series of data runs was taken by 
both adjusting heater power and auxiliary condenser coolant flow to obtain the same 
pressure conditions as above with a constant vapor velocity of ~1 m/s. 



49 



20000 




O 

CD 



O 

in 



o 



o 

CO 



o 

c\j 






CO 



I 




o 

o 

o 

CD 



O 

O 

o 

C\J 



o 

o 

o 

00 



o 

o 

o 



( >1 * Sviu / M ) OM 



Figure 5. Effect of Pressure and Vapor Velocity on the Steam 

Heat-Transfer Coefficient (Smooth Instrumented Tube, 
Heatex Insert) 
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Figure 6 shows how increased pressure, with vapor velocity held constant, gives 
a larger increase in the outside heat-transfer coefficient over the predicted Nusselt 
value. This discrepancy is also thought to be caused by the viscosity effect since the 
Nusselt treatment seemingly takes into account the other possible causes. Also 
depicted more clearly is the dramatic enhancement due to vapor shear effects for 
increased vapor velocity at constant pressure. 

As mentioned previously, the use of inserts allows more effective mixing of 
coolant and facilitates greater accuracy in the calculation of the outside heat-transfer 
coefficient. Figure 7 shows the outside heat-transfer coefficient determined from the 
use of the two inserts together with no insert. The results for the wire wrap and 
Heatex inserts are closely grouped, whereas the no insert data shows much greater 
scatter. This seems to indicate that the use of an effective insert does indeed 
enhance accuracy for the instrumented tube data. 

The mean temperature difference across the condensate film for the no insert 
case has a significantly lower value than either of the insert cases; with the outside 
and tube wall conditions unchanged, this indicates a higher inside resistance for the 
no insert case. To illustrate this point a mid-range instrumented tube data point at 
the same conditions for each insert case is shown in Table 3 (Pjgt—lOl kPa, 
Vvapor~ 11 m/s, V^,3„,~2.75 m/s). It shows that both inserts have a comparable effect 
(slightly better with Heatex) and roughly provide a factor of two enhancement over 
the no insert case in the inside heat transfer coefficient. Due to the increased inside 
resistance for the no insert case, the heat flux shows a decrease of about 15%. 



51 



14000 



O 

CD 





03 

CL 




CT3 

Q- 
















o 




CO 




03 






CVJ 


^ 


CL 


— 


03 




CC3 

CL 






Q- 




O 


o 

o 

x: 


CO 

CVJ 


o 

CD 

x: 


CO 

CVJ 




h- 


•* 


H 




in 




in 




CO 












E 


o 


E 


o5 


E 


f- 




q 


CO 

in 


q 






T“ 


3 


CO 




z 




z 




s 


1 

i 


< 




# 



# 



# 







o 

CVJ 



"T 

o 

o 

o 

CO 




( >t , 3vUi / M ) oq 



o 

o 

o 



Figure 6. Comparison of Experimental Results with Nusselt Theory 
for Varying Pressures and Vapor Velocities (Smooth 
Instrumented tube, Heatex Insert) 
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Figure 7. Comparison of Steam Heat-Transfer Coefficients at Atmospheric 

Conditions (Pj 3 ,— 101 kPa, m/s) for Three Insert Conditions 

(Smooth Instrumented Tube) 
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Table 3. PERFORMANCE COMPARISON OF TUBE INSERTS 



AT ATMOSPHERIC CONDITIONS 





Heatex Insert 


Wire Wrap Insert 


No Insert 


h, (kW/m^ K) 


30.85 


27.63 


15.31 


q (kW/m^) 


481.4 


477.7 


409.6 


h„ (kW/m' K) 


9.86 


9.82 


10.75 


AT(K) 


48.82 


48.65 


38.10 



Figure 8 shows the wall temperature profiles for the data points listed in Table 
3 along with a vacuum data set (Ps,,~28 kPa, m/s, V^^^n,~2.75 m/s) for 

comparison, where 0 degrees is at the top dead center position of the tube. The 
shape of the temperature profiles shows the effect of condensate film thickening 
toward the bottom of the tube. The higher resistance through a thicker condensate 
film results in a lower tube wall temperature toward the bottom of the tube as 
shown. 

As expected, the wire wrap insert and Heatex insert temperature profiles at 
atmospheric pressure are very similar. The no insert profile has a shape similar to 
the wire wrap and Heatex profiles, but the mean wall temperature is about 11 K 
higher and shows the effect of the increased inside resistance to heat transfer. The 
vacuum run temperature profile shows the effect of a lower temperature gradient 
between the steam and coolant; the lower heat-transfer potential results in lower 
heat fluxes and a flatter temperature profile. 

Figure 9 shows the comparison of the instrumented data with the predictions 
of Nusselt, Fujii, and Shekriladze and Gomelauri covered in section II B. The data 
depicted ranges from = 28 kPa and = 3.5 m/s to = 101 kPa and 
= 1.1 m/s. The data seems to follow the Shekriladze and Gomelauri prediction the 
closest, but is also very near the Fujii prediction. 
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Figure 8. Horizontal Tube Wall Temperature Profiles (Smooth 
Instrumented Tube, m/s) 
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Figure 9. Comparison of Experimental Results with the Predictions 
of Nusselt, Fujii, and Shekriladze-Gomelauri. 
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C. INSIDE HEAT-TRANSFER CORRELATIONS FROM INSTRUMENTED 



TUBE RESULTS 

The instrumented tube provides a method of determining an inside heat- 
transfer coefficient for each insert from direct wall temperature measurements. 
Wanniarachchi et al [Ref. 8], in a similar effort to resolve the differences between 
the Sieder-Tate correlation and his experimental results, developed a correlation 
based on the Sieder-Tate expression. Using an intercept form, with data taken on 
the same apparatus as used in this study, and a least-squares fitted leading 
coefficient, the correlation took the following form: 



Nu=0.0(ARe^^Pr^^ 



\0.14 



+ 26.4 



( 6 . 1 ) 



However, Rouk [Ref. 4], using an optimization technique, showed that the value of 
the intercept had little effect on the results, and it was the accurate determination 
of the Reynolds exponent that was more critical. 

With a view to finding the appropriate Reynolds exponent and leading 
coefficient for equation (5.23), the data for no insert, wire wrap, and Heatex insert 
were plotted as In Re versus In (Nu/Pr'^(^j/)ix,^)® ''') as explained in section V C. The 
plotted data are shown in Figures 10, 11 and 12 for no insert, wire wrap insert, and 
Heatex insert respectively. A line of best fit (typically with a regression coefficient 
of 0.99) was used to obtain the value of the intercept. Figure 13 shows the plot of 
all three cases on the same graph for comparison. Note that the insert data is 
closely grouped on the upper regions of the graph when compared to the no insert 
data. The increase in the Nusselt number again indicates more efficient inside heat 
transfer for the insert vice the no insert case. The Reynolds exponent, or slope. 
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Figure 10. Log-Log Plot of Re versus for No 

Insert (Smooth Instrumented Tube) 
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In Re 




Figure 11. 



Log-Lx)g Plot of Re versus for Wire Wrap 

Insert (Smooth Instrumented Tube) 
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Figure 12. Log-Log Plot of Re versus for Heatex 

Insert (Smooth Instrumented Tube) 
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10.0 10.5 11.0 



differs from the Sieder-Tate-type equation of 0.8 in all three cases. The following 
derived correlations apply specifically only to the medium tube, but should be 
applicable to any tube with the same inside diameter. 

The no insert inside heat-transfer correlation had the form: 



Equation (6.2) was used to reprocess some current smooth tube runs, with no insert, 
as well as those of Guttendorf [Ref. 3] and Van Petten [Ref. 2] to check their values 
for the smooth tube Nusselt coefficient a. 

The wire wrap insert inside heat- transfer correlation had the form: 



Equation (6.3) was used to reprocess previous data taken on the same apparatus by 
Guttendorf and Van Petten, who used the wire wrap insert for their finned and 
smooth tube experiments, with a view to checking their results with this new 
correlation. Table 4 shows a comparison between the new wire wrap insert 
correlation (eq. 6.3), and Wanniarachchi’s correlation (eq. 6.1) both with and without 
the intercept value included. This comparison shows that equation 6.1 is actually 
more accurate without the intercept value included when compared to the results of 
the new wire wrap insert correlation. 



ATu =0.013 — 




( 6 . 2 ) 



Nu= 0.052 — 



u \® 



(6.3) 
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Table 4. COMPARISON OF EQUATIONS (6.1) AND (6.3) 



FOR (held cons 


tant for comparison) 


Reynolds Number 


20,000 


30,000 


40,000 


50,000 


Nusselt Number 

Nu= 0.064 Re® «Pr‘'^(At7AtJ° ‘'-h26.4 


273.7 


368.4 


456.9 


541.0 


Nusselt Number 
Nu=0.064 


247.3 


342.0 


430.5 


514.6 


Nusselt Number 
Nu=0.052 


244.9 


341.5 


432.3 


519.1 



The Heatex insert inside heat-transfer correlation had the form: 



Nu =0.22 Re^^Pr^'^ 



/ \0,14 



(6.4) 



Equation (6.4) was used to reprocess all Heatex data. 

Memory [Ref. 12] conducted condensation experiments on a different apparatus 
with a smooth instrumented tube. He also determined the Reynolds number 
exponent for no insert, wire wrap insert (made locaUy and somewhat different from 
the one used in this study), and Heatex insert. The exponents he obtained are 
reported in Table 5. TTie Heatex insert exponent of 0.68 compares very favorably 
with the value of 0.69 found in this study. The wire wrap exponent of 0.73 was weU 
below the value of 0.82, most likely due to differences in the insert. The no insert 
case gave a value of 0.85 and compared well to the value of 0.89 found in this study. 
The difference in the no insert exponent is most likely due to the difference in the 
tube entrance region; Memory used a long run of straight pipe for the tube entrance, 
whereas the present work had a sharp bend just prior to the condenser tube. 

Rouk [Ref. 4] used the instrumented smooth tube data of Georgiadis [Ref. 5] 
to the find the appropriate value of the Reynolds number exponent for the no insert 
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and wire wrap cases. His results compare quite well with this study and are also 
given in Table 5. 

ANL (Argonne National Laboratory) [Ref. 21] conducted an assessment of 
heat-transfer correlations for turbulent pipe flow with water to determine the best 
correlation (s) on which to base their design of Ocean Thermal Energy Conversion 
(OTEC) heat exchangers. ANL used two shell-and-tube heat exchangers, with no 
inserts, for analysis and reported the following; 

1. The Dittus-Boelter (eq 2.7) and Sieder-Tate (eq 2.9) correlations under- 
predicted the data by 5% to 15% and were considered too conservative for 
design. 

2. Overall, the "best" correlations were found to be Petukhov-Popov (eq 6.5) and 
Sleicher-Rouse (eq 6.6), both of which showed excellent agreement (± 5%) 
with the experimental data (at Pr=6.0 and Pr=11.6). 



(e/8)jfePr 

(valid for 0.5 < Pr < 2000 and 10'' < Re < 5x10*^) 

where; 

€ = (1.82 logio Re - 1.64)-^ 

K, = \ + 3Ae 

K 2 = 11.7 + 1.8 Pr 

N« = 5+0.015 Re/Pr^ 

(vahd for 0.1 < Pr < 10^ and 10' < Re < lO"*) 

where; 

a = 0.88 - 0.24/(4+PrJ 

b = 1/3 + 0.5e-^-^'^ 



(6.5) 



( 6 . 6 ) 
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3. The most accurate correlations (i.e. Petukhov-Popov and Sleicher-Rouse) 
seem to yield effective Reynolds exponents in the neighborhood of 0.85 
(uncertainty range: m = 0.82 to 0.88). 

4. The potential sources of uncertainty in the Wilson procedure included 
waterside flow maldistribution, entrance effects, experimental error in U , and 
the uncertainty in the Reynolds number exponent. Of these, they concluded 
that the uncertainty in the Reynolds number exponent was, by far, the most 
significant. In fact, the results of the Wilson procedure were found to be 
highly sensitive to the value of the Reynolds number exponent. 



Table 5. COMPARISON OF REYNOLDS NUMBER EXPONENTS FOR 



SIEDER-TATE-TYPE CORRELATIONS 





Experimental 

Data 


Rouk 


Memory 


ANL 


No Insert 


0.89 


0.90 


0.85 


0.85 


Wire Wrap 


0.82 


0.78 


0.73 


— 


Heatex 


0.69 


— 


0.68 


— 



D. ANALYSIS OF SMOOTH TUBE RESULTS 

When using the modified Wilson plot technique to reprocess data files, the 
solution option can be specified to use either the stored value of the Sieder-Tate 
coefficient (for direct computation of h^) or let the coefficient value "float", which 
allows the program to calculate its own value of the coefficient. In order to 
determine which method was most accurate, the instrumented data files were 
reprocessed using each method and then compared with the values of the heat- 
transfer coefficient which were obtained by direct measurement of the tube wall 
temperature. A high, medium, and low coolant flow rate was chosen from each run 
to facilitate the comparison. The results were tabulated and are shown in Appendix 
D; it can be seen that the fixed coefficient method yielded the more accurate results 
at least 75% of the time. The mean error of the fixed method was ±2.0%, and that 
of the floating method was ±5.4%. The error for the lowest coolant flow rate (Re 



65 



<20,000) was noticeably higher than for higher coolant flow rates for both methods. 
The choice of using the fixed coefficient method represents a departure from the 
practice of previous researchers on this apparatus who exclusively used the floating 
coefficient method. 

Prior to the instrumented tube runs reported in Figure 5 (from which the new 
correlations were empirically derived) a series of data runs were made using a plain 
smooth uninstrumented tube (S02) of the same dimensions using the Heatex insert. 
The plain smooth tube data was then reprocessed using the new correlation for the 
Heatex insert and the results are plotted in Figure 14. These data sets were taken 
at the same conditions as the instrumented data of Figure 5 except for the set at the 
highest vapor velocity of 6.2 m/s vice 3.5 m/s for the instrumented tube. Similar 
effects of vapor shear and vapor pressure, as mentioned previously for Figure 5, are 
clearly seen, and again illustrate the vapor shear effect on the outside heat-transfer 
coefficient. 

With the exception of the two data runs at high vapor velocity, the data from 
Figures 5 and 14 are shown together in Figure 15. The close agreement of the 
reprocessed plain smooth tube data with the instrumented smooth tube data allows 
a high degree of confidence in the accuracy of the new correlations and the choice 
of the fixed coefficient method. 

To provide a baseline from which to evaluate finned tube performance it was 
necessary to obtain the smooth tube Nusselt coefficient, a, for the specific conditions 
under which the comparison was to be made. The condition chosen was atmospheric 
pressure (101 kPa) and a vapor velocity of —1 m/s. 

For 8 complete sets of data the average value of a, using the fixed method, was 
found to be 0.876. The average value for each data set was found by taking the 
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Figure 14. Effect of Pressure and Vapor Velocity on the Steam 
Heat-Transfer Coefficient (Non-instrumented Smooth 
Tube, Heatex Insert) 
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Figure 15. Comparison of Instrumented Smooth Tube Results with Non- 
Instrumented Smooth Tube Data After Reprocessing with the 
New Heatex Insert Inside Heat-Transfer Correlation. 
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measured value of h^, for each data point, dividing it by the Nusselt theory prediction 
of ho, and then multiplying by 0.728 (the Nusselt coefficient); the average of all the 
data points in the set was then taken. Interestingly the Wilson plot floating 
coefficient method gave a value of 0.835, somewhat lower than the average value. 
Originally it was thought that this discrepancy might be due to "outlier" data points 
(high or low coolant velocity) in each set. However, removing the highest or lowest 
coolant flow rates within a set had little or no effect on the Wilson plot result. The 
reason for the discrepancy is still not known and merits future study. 

The value of a was calculated for several other flow conditions; these are 
shown in Table 6. The trend of the readings, like that of Figure 6, shows that vapor 
velocity has a much greater effect than pressure on the value of a, as expected. 



Table 6. SMOOTH TUBE a SUMMARY; EFFECT OF PRESSURE AND 
VAPOR VELOCITY 



File Name 


P(kPa) 


Vapor Velocity 

(m/s) 


a 


FIMAVSHl 


28 


3.5 


1.015 


FIMAVSH2 


41 


2.5 


.985 


FIMAVSH3 


68 


1.5 


.930 


FIMAVSH4 


101 


1.1 


.866 


FIMAVSH7 


69 


1.0 


.836 


FIMAVSH6 


41 


1.0 


.818 


FIMAVSH5 


28 


1.0 


.786 



E. ANALYSIS OF nNNED TUBE RESULTS 

With an accurate value of a, and the newly determined inside heat-transfer 
correlations, the medium family finned tube data of Van Petten [Ref. 2] was 
evaluated. Figure 16 shows the data Van Petten reported in his thesis; it also shows 
his data after being reprocessed using the new wire wrap insert correlation with 
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Figure 16. Comparison of the Steam Heat-Transfer Enhancement Data of Van 
Petten, for the Medium Finned Tube Family, Using the Modified 
Wilson Plot and New Wire Wrap Insert Inside Correlation. 



70 



fin spacing ( mm 



both the fixed coefficient and Wilson plot floating coefficient methods. Since Van 
Petten used the Wilson plot method, it is not surprising that the original thesis data 
and new Wilson plot floating coefficient data are comparable since the Reynolds 
exponent only varied from 0.8 to 0.82. The fixed coefficient method enhancement 
is substantially higher than the Wilson plot results. 

Since the assertion is that the fixed coefficient method is more accurate than 
the Wilson plot method, then the conclusion must be that the enhancement for this 
set of finned tubes is actually higher than previously reported. 

During this study, limited medium finned tube experiments were conducted for 
purposes of comparison. Figure 17 shows the comparison between this data and the 
newly reprocessed data of Van Petten (using the fixed coefficient method) and shows 
reasonable agreement. To more clearly illustrate this point, the data taken on the 
2.0 mm fin spacing tube has been given in more detail in Figure 18. Excellent 
agreement is seen between the experimental results of this study and that of Van 
Petten using the known inside heat-transfer correlation with the fixed coefficient 
method. Again, as shown in Figure 16, the Wilson plot prediction is significantly 
below the fixed coefficient results. 
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Figure 17. Comparison of the Steam Heat-Transfer Enhancement Data 

of Van Petten and Swensen for the Medium Finned Tube Family. 
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Figure 18. Comparison of the Steam Heat Transfer Data of Van Petten and 
Swensen for the 2.0 mm Fin Spacing Medium Tube. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 



1. The inside heat-transfer correlation is highly sensitive to the Reynolds number 
exponent. 

2. Each insert condition must be analyzed separately to determine the 
appropriate "form" of the inside heat-transfer correlation. 

3. Calculations based on a known inside heat-transfer correlation are more 
accurate than modified Wilson plot results. 

4. Armed with accurate inside heat-transfer correlations, previous data may be 
reprocessed to give more accurate results. 

5. The source of contamination in the test apparatus, which has caused a 
dropwise condensation problem for a number of years, is most probably due 
to a contaminated distilled water source. 



B. RECOMMENDATIONS 



1. Reprocess all previous medium and large diameter finned tube data using the 
fixed coefficient method to obtain more accurate results. 

2. Continue with construction of smooth instrumented tubes of different 
diameters (i.e. small, medium, and large) to confirm the medium tube results 
and develop correlations specifically for the smaU and large diameter tubes. 

3. Construct one representative instrumented finned tube to test the validity of 
applying instrumented smooth tube results to finned tube data. 

4. Test representative water samples that have been collected from both the 
distiller and boiler to confirm the presence of impurities and validate their 
origin. 

5. Replace current distiller with a deionized pure water source (either 
commercial purchase of water or new distilhng apparatus). 
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APPENDIX A. PHYSICAL AND THERMODYNAMIC PROPERTIES OF 
WATER 

The physical and thermodynamic properties of water were based on the 
following equations: 

Heat Capacity (J/kg-K): 

Cp = 4.2 1 1 - r * [2.268X lO'^ - T * (4A24x lO'^ + 2.7 14x 10*’ * 7)] ( A. 1 ) 

where: T = temperature (celsius) 

Dynamic viscosity, (kg/m s): 



p=(2.4;cl0-5)*10f’^’»/(^^*”.i5)] 



(A.2) 



Thermal conductivity (W/m-K): 



-0.9225+x *(2.8395 -x *(1.8007 -x *(0.5258 -0.0734 *x))) 



(A.3) 



where: x = (T4-273.15)/273.15 



Density (kg/m^): 



p =999.5295 +T* (0.0127-7’* (5.4825x70*^-7’* 1.2341x70*5)) 



(A.4) 
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Latent heat of vaporization (J/kg): 



=2477200 -2450 * (T-10) 



(A.5) 



Fluid Enthalpy (J/kg): 



h^= T* (4.2038 - T* (5.8813;ci0 -^-7* 4.5516x70®)) 



(A.6) 



Saturation pressure (Pa): 



7^=22120000 

sat 



(A.7) 



where: 



Pr = e®' 

Br = SUM /[T, (1+ 4.1671 (1-T,) + 20.9751 (1-T,)^) - (1-T,)/ 

((1.0 X 10^(1-T,)^ + 6)] 

SUM = (-7.6912)(1-T,) - 26.0802 (1-T,)^ - 168.1707 (1-T,)^ + 64.2329 (1-T,)' - 
118.9646 (1-T,)' 

T, = (T+ 273. 15)/ 647.3 
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APPENDIX B. SYSTEM CALIBRATIONS AND CORRECTIONS 



B.l Thermocouple and Quartz Thermometer Calibration 

Several different thermocouples and the quartz crystal thermometer (HP 
2804A, Ser. No. 2244AD1192) were cahbrated against a platinum resistance probe 
in a mixed isothermal ethylene glycol bath with a Rosemont Galvanometer model 
920A commutating bridge (Ser. no. 013494) from 23 to 26 September, 1991. 

NOMENCLATURE 

T1 Quartz Thermometer Measuring Probe T1 (Ser. No. 2120A-00707, with dial 

setting 481 ) 

T2 Quartz Thermometer Measuring Probe T2 (Ser. No. 2120A-60459 with dial 

setting 510) 

T1-T2 The T1-T2 reading on the quartz thermometer 

T-55 Large diameter metal sheath thermocouple (diameter = 0.040") 

T-56 Small diameter metal sheath thermocouple (diameter = 0.020") 

T-57 Old Thermocouple (taken from rig during disassembly set to HP 3497A 
internal zero); (Type: Omega, TT-T-30, Lot# HCP093HC0306) 

T-58 New 1 Thermocouple (set to HP 3497 A internal zero); (Type: Omega, TT- 
T-30 SLE, Lot# OCP1453PTCC01473P) 

T-59 New 2 Thermocouple (referenced to ice bath zero) 

T-60 10 Junction Thermopile #1 (referenced to ice bath zero) 

T-61 10 Junction Thermopile #2 (referenced to ice bath zero) 

T-62 10 Junction Thermopile #3 (referenced to ice bath zero) 

The cahbration was performed by taking the bath temperature up from 290 K 
to 393 K by 5 K increments then back down to check for hysteresis; no hysteresis was 
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observed. The data was fitted to a fifth-order polynomial (regression coefficient = 
1.000 for each polynomial fit) in each case with the following results: 



T55 = 273.15+(2.5943c-2)K-(7.2671-7)T2 
+ (3.2941 -11)T^-(9.7119-16)K^ 

+ (9.7121-20)1^ 

T56 = 273.15+(2.5878c-2)K-(5.9853e-7)K2 

- (3.1242e-ll)K^+(1.3275e-14)K'‘ (B-2) 

- (L0188c-18)K^ 

T57 = 273.15+(2.5923e-2)K-(7.3933e-7)F2 

+ (2.8625c-ll)r^+(1.9717€-15)K^ (®-^) 

- (2.2486e-19)K^ 

T58 = 273.15 +(2.5931c-2)K-(7.5323c- 7)F2 
+ (4.0567c-ll)K^-(1.2791c-15)K^ 

+ (6.4402e-20)K^ 

T59 = 273.15+(2.5471c-2)K-(3.7621c-7)K^ 

- (1.0105c-10)T3+(2.3928e-14)K^ (®-^) 

- (1.6440e-10)l^ 

T60 = 273.17+(2.5571c-2)K-(1.9980e-7)T2 

- (1.6385e-10)K^+(2.6164e-14)F^ (B-6) 

- (1.0295^-18)1^ 

T61 = 273.15+(2.6119e-2)F-(9.0449e-7)T^ 

+ (1.1214e-10)F^-(1.5623e-14)F^ (®-2) 

+ (1.0646e-18)l^ 



78 



(B.8) 



T62 = 273.15+(2.5996c-2)F-(7.4405e-7)K2 
+ (2.4733<r-ll)F3+(3.3236e-15)F‘‘ 

- (3.7460e-19)F’ 

For T55, T56, T57, T59\ V= Voltage in microvolts (eg. 0.010 volts = 10,000 
microvolts). 

For T60, T61, T62; V= Voltage in microvolts/10 (eg. 100,000 microvolts ^ 10 

= 10 , 000 ). 

The HP 2804A Quartz Crystal Thermometer was also cahbrated on 26 
September, 1991 with the results summarized in Table B.l. 



Table B.l. QUARTZ THERMOMETER CALIBRATION DATA 



Reference 


T1 


Error 


T2 


Error 


Temperature 
(deg C) 


(deg C) 


(deg C) 


(deg C) 


(dec C) 


17.979 


18.003 


0.024 


18.001 


0.022 


18.977 


18.997 


0.020 


18.996 


0.021 


19.984 


20.001 


0.017 


20.002 


0.018 


20.980 


21.007 


0.027 


21.008 


0.028 


22.936 


22.959 


0.023 


22.961 


0.025 


25.014 


25.032 


0.018 


25.037 


0.023 


30.140 


30.153 


0.013 


30.164 


0.024 


34.977 


34.983 


0.006 


34.998 


0.021 



The mean error for the Quartz Thermometer as shown in the Table B.l is 
— +0.02 deg C. However, it is well noted that T1-T2, the critical measurement, 
indicates an apparent error of less than +0.005 deg C when the temperature reading 
is near or below 25 deg C; 77 and T2 tracking well together lowers the error 
estimate. 
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B.2 Flow Meter Calibration 

The flow meter calibration for coolant flow through the single horizontal tube 
was completed on Oct 28, 1991 using a stop watch, portable tank, and a Toledo 
model 31-0851 IV, Se. No. 1326 scale with 1/10 pound graduations. The following 
relation was obtained via linear regression: 



m = 6.7409F + 13.027 



B.9 



where: 

m = mass flow rate (grams per second) 

F = flow meter reading (eg. 10%=*>10) 

The applicable range of the calibration was 10% to 95%. The water 
temperature on the day of the calibration was 17.5°C (290.6 K) and water density was 
998.5 kg/m^. The data is shown in Figure B.l 

B.3 Mass Flow Rate Correction 

The inlet water temperature from the coohng water sump varies anywhere from 
15°C to 25°C depending on environmental conditions. To account for these 
temperature variations the following function was used to calculate a correction 
factor for viscosity variation with temperature [Ref. 5]. 

C/=1.0365 -(1.9644£-3)7In +(5.2500e-6)71rt2 

where: 

Cf = mass flow rate correction factor 
Tm = inlet temperature (celsius) 
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Figure B.l Horizontal Tube Coolant Flowmeter Calibration Chart 
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Flow % 



The value of Cf for the flow meter cahbration (Tm = 17.5°C) was 1.0037. 
Therefore, the actual mass flow rate was calculated using the following equation: 






Cf 



1.0037 



(B.ll) 



where: 

m^c, = corrected mass flow rate (grams per second) 
mcaic= computed mass flow rate (eq. B.9) (grams per second) 



B.4 Pressure Transducer Calibration 

Three methods of pressure measurement were available on the apparatus: 

1. Direct pressure reading off the Heise solid front - CM-104119 pressure gauge, 
(range 0-15 psia). 

2. Converted voltage readings from the Setra, model 204, Ser. no. 63982 
pressure transducer (range 0-14.7 psia; 0-5 volts; 5V—0 psia). 

3. Steam saturation temperature measurement with the apparatus producing 
steam at steady state. The steam saturation temperature/pressure relation 
was utilized via standard steam tables. 

The pressure transducer was cahbrated versus the vapor temperature probe 
reading on 12 December 1991. Equation B.12 gives the desired relationship. The 
data is shown in Figure B.2. 



P = -2.9360K+ 14.7827 



(B.12) 



where: 

P = pressure (psia) 

V = pressure transducer voltage reading (volts) 



B.5 FRICTION TEMPERATURE CORRECTION 

As coolant flows through the tube there is a bulk temperature rise due to 
frictional heating, which is highly dependent on fluid velocity. Although small, this 
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Figure B2 Pressure Transducer Calibration Chart 
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Transducer Voltage (volts) 



can have a significant effect on the calculated overall heat-transfer coefficient. 
Measurements were made for no insert, Heatex, wire wrap and twisted tape inserts 
as shown in Figure B.3, on 5 December 1991. Each data set was curve fitted to a 
third order polynomial which is depicted in Table B.2. Each respective polynomial 
corrects the temperature rise measurement for the heating due to the particular type 
of insert used. 



Table B.2 FRICTION TEMPERATURE RISE POLYNOMIALS 



Insert Tvoe 


Polynomial 


None 


Trise = -1.960x10 ^V^-l-9.349xl0''V2-Kl.749xl0“V-2.728xl0'' 


Wire Wrap 


Trise = 8.160x10 "'V'-H.4512xl0-^V2-i-2.745xl0-^V-3.991xl0“ 


Heatex 


Trise = 8.160xl0-^V^-t-1.080xl0-^V2-t-1.232xl0V 4-8.570x10-^ 


Twisted Tape 


Trise = 4.070x10 ^V^-f4.451xlO-''V^-»- 1.71 1x10 ^V-6.440xl0^ 


where: Trise 
V = 


= temperature rise (K) 
fluid velocity (m/s) 
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Figure B.3 Friction Temperature Rise Curves for Heatex Insert, Wire Wrap 
Insert, Twisted Tape Insert, and No Insert. 
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Flow Velocity (m/s) 



APPENDIX C. SYSTEM INTEGRITY / LEAK TESTING 

As mentioned in section III C, the material for the auxiliary condenser 
penetration plates was changed from aluminum to stainless steel on 24 January 1992. 
The stainless steel screw thread connectors in the aluminum side plate had loosened 
to the extent that leakage could be detected. The cause of this loosening was due 
to thermal cycling of the apparatus and the differential contraction/expansion of the 
aluminum/stainless steel combination. 

An initial leak test was conducted from 20 December 1991 through 2 January 
1992; the results are shown in Figure C.l. The initial mean leak rate was 3.4 mmHg 
per day. 

Subsequent to the structural modification noted above, another leak test was 
conducted 6-19 February 1992; the results are shown in Figure C.2. The mean leak 
rate was found to be 1.7 mmHg per day, an noticeable improvement. In general, a 
leak rate of 2 mmHg per day is considered acceptable. 
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Figure C.l Apparatus Leak Test I 
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Figure C.2 Apparatus Leak Test II 
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Time (days) 



APPENDIX D. COMPARISON OF FIXED Q vs FLOATING C, SOLUTION 
METHODS FOR MODIFIED WILSON PLOT DATA REPROCESSING. 



As related in section VI D the "fixed" C, method and "floating" C, method were 
evaluated against the original instrumented tube results in order to choose the most 
accurate method of reprocessing non-instrumented data. The results of the fixed 
coefficient versus modified Wilson plot floating coefficient method comparison are 
shown in Table D. The fixed coefficient method was determined to be more accurate 
than the floating coefficient method 75% of the time. For the three coolant flow 
rates considered for each run, the overall fixed coefficient method mean error was 
±2.0%, and the floating coefficient method mean error was ±5.4%. 

For Re <20,000 the error was noticeably higher, which tends to support the 
assertion of a number of researchers that data with Reynolds numbers below 20,000 
should not be used, since the flow may not be fully turbulent. 
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Table D. COMPARISON OF FIXED C, vs FLOATING C, REPROCESSING 
METHODS 



File Name 


Flow 

Vel. 


Film 

AT 


Inst 

Tube 

Data 


Fixed 

Method 


% 

error 


Float 

Method 


% 

error 




V. 

(m/s) 


(fe 


(kW/m^-K) 


ho 




ho 




FIMAVSHl 


4.35 


29.33 


12.13 


11.93 


-1.7 


13.19 


8.7 




2.21 


22.81 


13.18 


12.82 


-2.6 


15.17 


15.1 




1.16 


18.27 


13.68 


13.16 


-3.8 


17.20 


25.7 












-2.7* 




+ 16.5 


FIMAVSH2 


4.34 


36.09 


11.10 


11.14 


0.4 


11.07 


0.3 




2.22 


30.39 


12.03 


11.87 


-1.3 


11.72 


-2.6 




1.16 


24.62 


12.66 


12.50 


-1.3 


12.26 


-3.1 












1 

o 

* 




-1.8 


FIMAVSH3 


4.35 


46.10 


10.31 


10.26 


-0.5 


10.62 


3.0 




2.23 


39.35 


10.89 


10.83 


-0.6 


11.44 


5.1 




1.16 


32.56 


11.42 


11.40 


-0.2 


12.43 


M 












* 

o 

1 




-t-5.6 


FIMAVSH4 


4.35 


55.20 


9.263 


9.146 


-1.3 


9.584 


3.5 




2.23 


47.31 


9.869 


9.756 


-1.1 


10.52 


6.6 




1.16 


39.45 


10.41 


10.39 


-0.2 


11.71 


12.5 












-0.9* 




-t-7.5 


FIMAVSH5 


4.35 


32.83 


8.618 


8.478 


-1.6 


9.018 


4.6 




2.23 


28.05 


9.589 


9.341 


-2.6 


10.36, 


8.0 




1.16 


22.99 


10.23 


9.865 


-3.6 
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File Name 
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File Name 


Flow 

Vel. 


Film 

AT 
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Method 


% 
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Method 


% 
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(m/s) 




**<>2 
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ho 
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2.0 
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-7.9 
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” 
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-6.2 


FIMASN6 
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0.2 
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11.77 


03 



+7.4 +4.9* 



92 



APPENDIX E. UNCERTAINTY ANALYSIS 



In the measurement of a physical quantity, there will always be a difference 
between the measured value of the quantity and the actual value. The magnitude 
of this difference depends on the accuracy of the measuring device calibration, 
operator experience, environmental effects, etc. Eventhough the error associated 
with a single measurement may be rather small, the error may grow to substantial 
proportions when combined with other measured quantities in a given calculation 
scheme. The best estimate of the difference between a calculated or measured 
quantity and the actual value of the quantity is known as the uncertainty. 

The uncertainty may be estimated by the method of Kline and McClintock 
[Ref. 24]. This states that for a quantity R, which is a function of several measured 
quantities (R=R(x,,X 2 ,X 3 ,...,x„)), the uncertainty in R is given by the following 
relation: 






IB. 

dx. 



/ 



w. 



dx. 



\2 



dx, • 



1/2 



(E.l) 



where: 

Wr = the uncertainty of the desired dependent variable 
Xi,X 2 ,Xj,...,x„ = the measured independent variables 
Wi,W 2 ,W 3 ,...,W„ = the uncertainties in the measured variables. 



Georgiadis [Ref. 5] gives a complete description of the uncertainty analysis for 
this experiment. The uncertainty analysis program written by Mitrou [Ref. 25] was 
used to calculate uncertainties. The uncertainty for runs using an insert ranged from 
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±2.2% to ±4.6%. The uncertainty for runs without an insert ranged from ±10.9% 
to ±17.1%. Sample outputs of the uncertainty evaluations are included in this 
appendix. 
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DATA FOR THE UNCERTAINTY ANALYSIS: 



File Name: 

Pressure Condition : 
Uapor Tenperaturc 



Water 

LU+ -r 



Flow Rati 
Ucloc i ty 



FNMAUSH4 
Atnospher 1C 



v‘a 

Heat Flux 
Tulje“netal 
Sieder-Tate constant 
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go ^ 779 
50.00 
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0.2200 



(Deg C) 
( n / 5 ) 

( W/n"2 ) 

( W/n.K > 



UNCERTAINTY ANALYSIS: 



VARIABLE 



PERCENT UNCERTAINTY 



Mass Flow Rate, Nd 1.27 
Reynolds Nur.ber , Re 1.4S 
Heat Flux , q 1.37 
Log-Mean-Ten Diff, LMTD .27 
Wall Resistance, Rw 4.40 
Overall H.T.C. , Uo 1.40 
Water-Side H.T.C. , Hi 1 .55 
Uapor-Side H.T.C. , Ho 2. 64 



DATA 
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97 



FOR THE UNCERT^^INTY F.HF.LYSIS^ 
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DATA FOR THE UNCERTAINTY ANALYSIS: 



File Nane: 


FNNAUSHG 






Pressure Condition: 
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Uapor Temperature 
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DATA FOR THE UNCERTAINTY ANALYSI 
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FOR THE UNCERTi^INTY ftNFiLYSIS: 



File Nane: 


FNNAF201 






Pressure Condition: 


At mespher i 


0(101 kPa ) 




Uapor Temperature 


- 


99.787 


< Deg ' 


Water Flow Rate (%) 


m 


80.00 




Water Ueiocity 


- 


A -a.-^ 

-T . ta/ 


( m / s ) 


Heat Flux 


as 


1 .4315+05 


{ IJ / irj 2 ) 


Tuhe-metal thermal cc 


nduc. 


395.0 


( W/m.K ) 


S iede, — Tate constant 




0.2200 





UHCERTftlNTY ANALYSIS: 



N^p T 



PERCENT UNCERTAINTY 



Ma«c Floui RsvC r 
Reynolds Nunher , 
U-.+ Flux 

I O'* ■PJ 

U.W^ IWII I 

Wal 1 Resistance , 

1 1 u T n I 

WVWILU 4 .X I ^ 

U'ater~S:de H.T.C. 

I u T r 

vi_ipwi Fit < *wi 



H: 



80 

14 

QC 

20 

40 

O'? 

%j I 

34 

9D 



RATA pnp Tuc U^iCERTAINTY ANALYSIS: 



F 1 i “ Nane : 
Pro; 

Umry, 



CMMAC9^ 1 



A 4 

. I i, . I 



; (101 K P a ) 
100.076 



/ n.r. r ^ 



U' a t 
Meat FI 



lla 1 + w 

VWJ-WV-^ vjr 



hemal con 



1 ^ IP^ r 4 I I 



> 1 eder-T a t ^ 



2 . 74 

0^041 P+0C 

5S5.0 

0.2200 



(U/n"2 ) 
n,\ 1/ \ 



I IMPCDT A T MTV 

WI«WU-I\IIIJ.I«I I 



ANALYS 1 5 : 



UAR I ABLE 



PERCENT U.NCERTftlNTY 



Nsss Flow Rate, Nd 
1 D. 

ivwjr<iwJ.Ui«,/ i<it_iiii_i^i y i>'w 

Uaa 4 - C 1 , .w 

Lcg-Nean“Ten Diff, LNTD 
Wall Resistance, Rw 
Overall H.T.C. , Uo 

u»4 , r4« u T r u; 

n.I.W.y 11 ^ 

Uapor-Side H.T.C. . Ho 



1.27 

I .52 
1 .35 
.14 

9 C"7 

A. « w • 

1 “7C 

I ■ 

1 .57 
2.65 
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n^TA pQp yup UNCERTAINTY ANALYS 


I S : 




Filz Nar.e: 

PrcsGLirc CoDditicn: 


FNNAF201 
At rnespher 


1C (101 kPa ) 




Uapcr T cnpera tijrc 


- 


100. 142 


(D 


Ulster Flew Rate (%) 


- 


20.00 




Water Velocity 


- 


1.15 


(rn/ 


Heat Flux 


- 


B . 405E+05 


( W / n 


T'jbe“'rietal thernai c 


endue . 


395.0 


( W/n 


Sicder-Tatc conatant 


- 


(71 7Gf7ir?k 

xt/ . 1 - ^ XL/ XL/ 





UNCERTAINTY ANALYSIS: 

u Ap T ADj_r 



PERCENT UNCERTAINTY 



Mb. 






c: 


1 ^... 




D 




1 I V/ 






' 


▲ W Uf 




' ' 




Dbx/ 






1 




M 


, , 


txn i 


1 X V- 7 


‘ ' 




X 


Ui ^ 


1 ^ (-4 


1 M 


Ubb 






IT 


1 . 








1 • w u» 






' 


i Ui 


J 




M 


1 


_ 


M 




^ ^ ^ 


T 




»VX 


l_ wy 








Cdl 1 1 


' 




* • 


I.Ib 1 






D 


£Sk C ^ 








xnr x_i X 


X 




* X 


W ^ -4. 




'' 




Ox/B 




B 


1 


1 H 




T 


1 


w V W 


' 


*-* 


■*■ 


• 


' 




U's t 








C X u 






u 






' 




X XJ 






' ' 


M ar^ 


rx 






Sid 


A 




u 


V 0>p 


U 














HAT A 




crnD 


TUtr 


U/ 1 1 1 








W 1 ■' 








C , 1 






N 










1 X X 








3i)e 


* 






D^b 






, , 


^ a. 


r 






1 X- 








1 w 






' ‘ 


1 i ar^ 












B 




V X-' H 










H 


— 




l.la 4- 


B 






C 1 ^ 


, , 




D 


xr» W V 








1 ^ w w 




1 \ 1 


1.1 a •» 


B 






1 J B 1 








vxr «_> L 








\/ W X 


w 








4- 




c 


lux 








1 1 w w 


X* 




' 










Tub 






IVX 


B 4- a 










C 




* ' 


w V x_i 








Sie 


u 


e 


r 


“Ta 


+ 


£ 





P a 4 A N Cj 



n , ^ f < MTH 
a D, , 

I w , I X 

' lln 



0 1 
1 3 
05 

ChQ 

W 

57 

05 

"7f7l 

t w 

CQ 



FSCNNASN 1 

A t Ti G = h w n c ( 10] K P 5 ) 
= I 00 , 1 55 

s Qf?i r?kf7i 



PtwGr"i*l CGr7ClwG 



35 

5495*^05 

(h 



{ n* 

X U/ V- 

( 1*^ / 5 

( U! / n '' 
(U'/rj. 



UNCERTAINTY ANALYSI: 



NAD T j^9{_r 




PERCENT 


UNCERTAINTY 


Mass Flew Rate 


, Md 




0 


. 90 


DBx<ir-k/-x1pJ. KIiiixaKb 

ixxxjr IIWxLl^ 

Ub a 4- r 1 , 

1 1WX_. L 1 X U/-V , H 


r, Re 




1 

1 


1 ] 
r?iG 

• W 4b 


1 » »^-T B.VX n 

i_ wy 1 .ww. . . w. , 


iff LHTD 






. 45 


Wall Resistanc 


s, Ruj 




2 


. 57 


Over a 11 H . T . C . 


Uo 




1 


. 1 1 


U'atsr-Side H.T 


.'c. , Hi 




15 


.41 


Uspor-Side H.T 


.C. . Ho 




10 


.91 
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DATA FOR THE UNCERTAINTY ANALYSIS: 
File Nane: FSONNASN1 



Pressure Condition: 


At nospher i 


c ( 101 kPa ) 




Uapor Tenperature 


- 


99 . 842 


(D 


Water Flow Rate (%) 


- 


cfTi aa 




Water Uclocity 


- 


2 . 75 


( r,/ 


Heat Flux 




4 .0125+05 


( W/n 


Tube-netal thernai c 


onduc. 


385.0 


( W/n 


Sieder~Tate constant 


* 


0.0130 





ijKtrcpjAj^TY ANALYSIS: 
UARIABLE 



PERCENT UNCERTAINTY 



Mace CT1/-M.I Mr^ 

I XwiM tiLJ 

Reynolds Number, Re 
Heat Flux, q 

4_Wy l.W^.I . W. . iJil . , L_I.ILJ 



1 27 
1 .£0 
1 .39 
.52 



Wall Resist 


D... 

*-• 1 iC W , 1 MA, 


2 . £7 






1 1 U T 

W V W 1 !_* J. X II. 1 


.C. , Uo 


1 A 7 

1 . “T i. 






Wat er ~S i de 


H.T.C. , Hi 


IS. 44 






vi^pwi ^xu<w 


H.T.C. Ho 


17.07 






RATA pnp Tuc ij[vjrppTA|TMTY ANALYSIS: 






File Nans: 


FSONMASN 1 








Pressure Cc 


ndit^'^n: Atno**'^^‘**r 


, ^ / 1 R 1 L. D = ^ 

X W \ 1 W 1 1 . 1 l_l / 






U Wl 1 Wl ip w 


a 4 - . ip<> sr 


100.035 


( 


Deg C 


Water Flow 


R a t e ( % ) “ 


7 f 7 i aCA 

1 ^ ■ W V' 






Water U e 1 c c 


4 w = 

X V J 


3.91 


{ «r, 
V 11 


/s ) 


Heat Flux 


= 


4.39EE+05 


( w / 


r."2 ) 


J w 1 ‘ 1 W ^ X 


thernal conduc. 


385.0 


( w/ 


n.K ) 



Sieder-Tatt 



0.0130 



UNCERTAINTY ANALYSIS: 



I lADT ADI D 



PERCENT UNCERTAINTY 



Mass Flow Rate , Md 
Reynolds Nwnber , Re 
Heat Flux , q 
Lcg-'Mean-Ten Diff, LMTD 

Uall Resistance, Rw 
* 1 1 u T r 11^ 

W»WI LJJ.A , WW 



. + c , 



■ C. 



'v'spor-Side H.T.C. 



Hi 

Ho 



•5 

Lm 

1 

15 

1 o 



21 

08 

41 
G7 
17 

42 
42 
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APPENDIX F. INSTRUMENTED TUBE CONSTRUCTION 



The medium diameter instrumented tube of Poole [Ref. 6], with six wall 
thermocouples, was fabricated by welding together a copper tube in three pieces 
after imbedding capillary tubes at mid depth in the tube wall. Teflon thermocouples 
were inserted into the capUIary tubes for tube wall temperature measurement. 

An attempt was made during this study to construct instrumented tubes of 
small, medium, and large diameter in the following steps: 

1. Take thick-walled copper tube bar-stock and machine to specified inside 
diameter, outside diameter, and length. 

2. Cur four evenly spaced slots, of given depth and width to accommodate metal 
sheathed thermocouples, from halfway along the tube to the end. 

3. Solder metal sheathed thermocouples into groves. 

4. Cur copper strips from another tube to fit into the top of the slot. 

5. Clamp copper strips to the slots (using jubilee clips) and solder in place. 

6. Turn off excess copper from strips to original outside diameter. 

7. Send tubes to plating shop and plate the whole tube to a given thickness of 
plate. 

8. Return tube to NFS machine shop and machine back down to desired outside 
diameter. 

The process was completed through step 2 of the above procedure and the 
tubes are ready for step 3. Several of the latter steps were attempted using a 
practice tube with monel wire placed in the slots on the tube vice thermocouples; this 
was to ensure the process was safe, and to prevent destruction of assets (small, 
medium, and large diameter tubes machined to specifications, and the associated 
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metal sheathed thermocouples). The following summary documents the steps taken 
and lessons learned. 



1. Completed tube fabrication through step 2 of the fabrication procedure 
(above). A practice tube was used from this point to continue the procedure. 

2. Cutting the copper strips proved difficult, so the monel (similar to the 
stainless steel sheath on the thermocouples) wires were silver soldered after 
being pinned in place. No lead was permitted in the solder since copper will 
not plate to lead, but will plate to silver. 

3. The practice tube was sent to a local contractor for electroplating. (Bay 
Custom Chrome, Marina, Ca.) 

4. Two different plating procedures were used on the tube: 

a. First, a cyanide bath treatment was used to electroplate the tube. 
However, this procedure severely scorched the surface of the tube, 
resulting in the return of the tube to the NFS machine shop for repair. 

b. After repair, the tube was again sent off and subsequently treated with 
an acid bath procedure, which resulted in an acceptable tube surface 
finish. 

5. The original thickness of electroplate was not thick enough for our 
application, so another acid bath treatment was performed to bring the 
surface thickness to the desired level. 

6. The final product was suitable for our intended application with the exception 
of the following items which need to be corrected: 

a. Existing voids (from non-uniform soldering) or flaws in the tube surface 
would not fill in to give a uniform thickness around the tube. 

b. The practice wires, tube end pieces, and tube interior were not 
protected properly during the electroplating process, resulting in 
unacceptable copper deposition on critical areas. 



It is recommended that the instrumented tube efforts be continued with 
particular attention to the following items: 



1. Silver soldering is an acceptable method of fixing the thermocouples in place, 
since the thermocouple melting temperature is ~1700°C (silver solder is 
applied at ~1100°C). 
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2. If the voids left in previous soldering efforts persist, then cut copper strips to 
place over the thermocouples to hold them m place. 

3. Ensure that critical areas on the tube are well protected during the 
electroplating process. 
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APPENDIX G. DRPINST PROGRAM LISTING 



The program DRPINST, which was used to collect all instrumented data, is 
listed in this appendix. 
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1 (TiOiCh I 

\ 00B ^ 

1 ChChD I 

t : 

1 00 Q ! 
1 (71 1 f7l I 
1 0 1 I ! 
1 f7i 1 r> I 

I I ^ 

I 0 ] 3 I 
\ a^ A \ 

I U/ I -T 

1015' 

I 0 1 E ' 

1017' 

1 (7i 1 D I 

I W/ I w 

1 r?> 1 Q I 

I w I ^ 

1020' 
1 (?»'? 1 I 

I ^ I 

1 t: I 

I w 

1024' 

1 025 I 
1 025 I 
1 02^^ ' 
1 029 ' 
1 02 9 ' 
‘ 030 ' 

1 03 1 I 

1 032 I 

1034' 

1 I 

•1 ay.c I 

1 037 ' 
1 039 ' 
10401 

1041' 
104^1 
1 CTiAC I 

1 ^4£ ' 

^ OiA’7 \ 

W T I 

1049' 
^049' 
I 050 ! 

1 fTiC 1 i 

I w ^ I 

1 0C0 » 

1 053 ‘ 

1 acA I 

^ 055 ‘ 

1 0CC I 



DRP I NST 

WRITTEN FOR INSTRUMENTED TUBES 
Dv MEMORY ( 1 0TH DECEMBER 1991 ) 



USED TO CO'^'^ECT DATft ON THE INSTRUMENTED 
TUBE FftBRICftTED BY POOLE (1983) FROM WHICH 
NEW INSIDE HEAT-TRftNFER CORRELf^TIONS 
EMPIRICALLY DETERMINED FOR NO INSERT 

up£.p JMCCDT P|MQ UCAJCV JMCCpT CASES 

K. S'WENSEN (FEB 1892 ). 



'WERE 
, WIRE 

D V 



PjLL INSTR'JMENT C A'.» I BRAT I ONS FOR FLO'W METER 
THERMOCO'JPLES THERMOPILE PRESS'JRE TRANSD'JCER 
FP I CT I ON TEMP RISE FOR NO INSERT 'WIRE 'WRAP INSERT 
HEATEX INSERT AND T'WISTED TAPE INSERT 'WERE 
PERFORMED BY SWENSEN ( SEP-DEC 1991 ). 

UA^IO'JS PROGRAMS 'WERE 'USED TO REPROCESS DATA 
CO't_',_ECTED ; THE APP'i_" CABLE FILES 'WITH THE APPROPIATE 

DCDDOrcrCC T KiC DOnrZDAMC ADCr a: 



FILES 



orDDorcrc 



TM 0 ppncpAM 



S'WEN 



[TK^I CTTMAIiCUl-'? 



DRP I NST2 



CUCTMCCM CTTMACNM C 
cue M c r M r T M A I J C U 3 _ c 

S'WENSEN NON“ I N^TR'JMENTED DATA 
(ALL FN .... FILES) 

G 'J T T E N D 0 R F D A T A 

MAM P^TJcm hat a 



npp T KJCT 1 

HDD T MCT 1 
HDDL C 

^ I I w 

npppc2 

DRPKS2 



FILE 



M AMC 



DEFINITIONS 



FOR SWENSEN’. 



LETTER /.NUMBER 
PCS I T I ON 



DEF I N I T ION 



1 ST 

2ND 

3RD 

4TH 



5TH 



CONDENSATION CONDITION 
D“DRCPWISE F*=FILMWISE 
T'JBE TYPE I “ I NSTRUMENTED 
Nj = MnM- j NSTR'JMENTED 
T'JBE TYPE M“MEDI'JM 
PRESSURE CONDITION 
A^ATMnpj-jpRJC (]0i kPA) 

'v'='JACUUM 

'WHEN A 2. U APPE.AR TOGETHER AS IN 
FIMAUSH1 IT REFERS TO A SERIES 
OF RUNS AT UARIOUS PRESSURES AND 
VAPOR 'JELOCITIES. 

TUBE TYPE S“SMOOTH 
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1 0C7 ! 

105SI 

1 \ 
I ^ t 

] 0E0 ! 

1 0 E n 
1002 ! 
1 0E3 ! 
10S4 ! 
I0G5! 
10GE i 
1 0E7 ' 
10E8! 
10S9! 
) 070 ! 
1 07 1 • 
1 072 ^ 

1 073 i 
1 074 I 
1 anc I 

I W I wl 

]07S ! 
1 077 ' 

1 a^Q I 

i w I w ' 
1070! 
1 OiQCTi I 

I W X/ 

1 0 Q 1 i 
1 OiQ'? 

t Cf’ W 4. 

1 

I V/ ^ ^ 

1 004 
1 0QC 
1 000 
1 007 
1 002 
1 

t w ^ w 

] 0S0 



F-FINNED 

GTH TYPE OF INSERT INSTALLED 

N-NONE 
Ul-UIRE DRAP 
H-HEATEX 

THE HEATEX INSERT WAS USED EXCLUSIVELY 
FOR FINNED TUBE DATA, SO POSITIONS 
G &• 7 REFER TO THE TUBE FIN SPACING 
05-0 . 5nm 
10-1. 0nr5 
15-1. 5nn 
20-2 . 0nn 

FINAL POSITION RUN NUMBER 

MEANING OF ALL FLAGS IN PROGRAM 

IFT: FLUID TYPE 

ISO: OPTION WITHIN PROGRAM 

IM: INPUT MODE 

TC0; FINNED OR SMOOTH 

INN.* INSERT TYPE 

lUIT.* LOOP NO. WITHIN PROGRAM 

TTnc- T(ic^cr nTAMCTCD 

TDQ; PRESSURE CONDITION 

Ton* niiTDiiT DcniiTDcn 

± V • \J\^ I t ^ I 

COM /Cc/ 0(7) 

COM /0c55/ T55(5) 

COM /Cc5E/ T5E(E) 

COM /0c57/ T57(5) 

COM /Cc59/ T59(5) 

COM /Fid/ Ift 

HTM T.,^C\ T,.WC\ 

QATA 0.10026091 25727. 943E9 “7G7245.S295 79025595.91 
HAjjA -0247426599 E.97E99E+11 -2 . EG 1 G2E+ 1 3 3.94079E+14 



1 1 
1 w 1 


READ 


C( ) 








1 092 


RATA 

1 1 I 1 1 


273. 15, 


2 . 5945E-2 


-7 ^ 


2E71E-7 ,3.2941E”1 1 


1 A 

1 


READ 


TCC / M \ 

( ^ ^ \ t 








1 094 


DATA 


9 7 •7. 1C 
1 ^ . 1 ^ , 


2 ^ c0726~2 


-5 . 


9S53E-7 ~3.1242E~1 


1 ^^Qc; 

1 \y 


READ 


T5E( * ) 








1 09E 


DATA 


'?'77. 1C 

1 wl • 1 ^ J 


7 CQ7’^C-'? 

U. , 


-7 ^ 


3933E~7 ,2 . 9E25E~ 1 1 


1 ^Q'? 

1 .j 1 


READ 


T57( * ) 








1 099 


DATA 


273.15, 


2 .5931 E“2 


_7 ^ 


C9T:7C_'7 4 ^CC^C-ll 


1099 


READ 


TCQ/ * S 

1 W V / 








1 1 00 


O c « r^ 


- . 1524 


1 T 

JL. f 1 JL. ^ ^ 


die 


neter of stainless 


1 1 1 
t i %J t 


Ax-PI^Ossp^Z 


/4 






1102 


L - . 1 


3335 


! Ccndsns 


ing 


length 


1 !33 


LI-.' 


0G0325 


1 Inlet c 


nd 


“fin length" 


1 104 


L 2 - . 


034925 


! Cutlet 


end 


"fin length" 


1 105 


PRINTER IS 1 








1 10S 


BEEP 











Lon 
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1 



V/ 1 


INP'JT 


*' c r 1 c p T 

t— l_ W 1 


FL 


1 ITH 

W X w 


( 0=CJPAM 


09 


PR I NT 


IICTMC " 

W vJ ^ 1 <( W 


4X , 


" "Sc 


lect optic 


rTlQ 


PRINT 


tlCTMC " 

w w A < 4 w 


BX 


.. .. ^ 


Take data 


1 0 


PRINT 


USING " 


C V 

w/t , 


" " 2 


Print raui 


1 1 


PRINT 


USING " 


C V 


" " 3 


Purge 



1 I 2 

1 13 
1 14 
1 1 7 

120 ' 
1 3 
1 25 

1 2 g 

1 3*? 

1 35 

1 39 
1 4 1 
144 
147 
1 50 

1 C 3 
1 cc 

I ^ o 

1 7 ) 

1 74 
1 77 

I C*:; 

1 OQ 
1 

1 QC 

1 QC' 

20 ^ 
204 
207 

2 1 0 
2 1 1 

1 -5: 



or re-process previous 
data (NOT C0NP'u.ETE ) 



TMpJJJ T.^ 

T r T c: ^ M 

XI X O I 

PRINTER IS 70] 

T j Qb»0 
ri C£iC 'l(?\Q 

W l_ l_ I I I \ I W w/ 

BEEP 

INPUT "ENTER MONTH DP.TE AND TIME ( MM : DD ^ HH ; MM t SS ) " Date$ 
OUTPUT 709; " TO" ; Dates 

Hi !TDI IT . '' T H “ 



T CT T , 1 TUCM 

BEEP 

INPUT "SKIP PRGE P.ND HIT ENTER" OK 



CMTCD '7^710 . n=» 4 . ^(C 
DDTMT " 






^0 tine 



. n«, 4 - 



DRP INST 



pcTMT IJCTKJ 0 " 1 0 X ""NOTE- Progran rane 
IF IjOb”! THEN 1199 
BEEP 

TMDMT " 9 HTrp TKipijT MQnc ( npTA 2-EXISTING FILE)" In 

BEEP 

TC 

BEEP 

INPUT "GIUE P NAME FOR THE FPU DPTr. FILE" D fileS 
Dp]MT ncTMc "icy nane ^ "" 14P";D fileS 

Dn^T n f ^ \ ic 



: T CM . 1 



TH n p , 1 ««t 



TMpijT "CMTcc GEOMETRY CODE i I'FINNED 0*PLPIN)" I f g 

ddtmtcd tc 1 



a MCMC / nrr Ai n t s " 

xj iiwi—i / 

1 TWISTED TPPE" 

'7 LIT DC 

4- WAi\u- yviviu 

-7 UCATCV" 



1 2 1 E 


BEEP 




1 p 1 q 


pp I MT 


" SE 


1 220 


PRINT 


" I 


1221 


PRINT 


" 


1 222 


P R I N T 




1 r p — 

• w 


PRINT 


" 


1 


T [^jpij J 


I nn 


1 225 


nijTpijT 


@F 


1 227 


T MpiJT 


" NO 


1 p 2 9 


IF Ifg 


-1 


1 229 


I NPUT 


" EN 


1 *^30 


eese 




1 9- 1 


c v*nB r?i 





TM I.IAI I '• T .4- 

J.l') yXl.UUV 
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1232 Fp-0 

1234 Fw*0 

1235 END IF 



1 237 

1 243 

1246 

1 249 
1 292 
1 255 
files 
1 259 
1251 

\ OCA 

I LJ -r 

1 267 

1 ■7'7(7l 

t ^ i 

\ 273 
I 276 
I 292 
I 295 
I 295 
I 349 
I 35 1 

1 aC a 

I ^ w/ -r 

1 357 

1 -XCHi 

I ^ LJ W 

1 aC*3: 

I Ur w 

1 3 CC 

1369 
1 37 c 

1 379 
1 459 

1 402 
1 403 

1465 

1 ACC 

I "T ^ LJ 

1 407 

1 472 
1 474 

1 477 
1 479 
1 479 

1 A Q{A 



OUTPUT @Filc;Inwt ,Fp,Fw,Fh 
ELSE 

IF Ijob-1 THEN 1255 
BEEP 

INPUT "GIUE THE NAME OF THE EXISTING DATA FILE " ,D_f 1 1 e$ 
PRINT USING “l 6 X,""Thi 5 analysis was performed for data 

IF Ijob-1 THEN 1270 

Nrun“ 1 9 

BEEP 

INPUT "ENTER NUMBER OF DATA SETS ( DEF- 1 9 ) " ,Nrun 

ASSIGN @File TO D fiieS 

ENTER @Fiie;Ifg,Inn 

ENTER @File;Inwt ,Fp,Fw,Fh 

END IF 

IF Ijob“l THEN 1543 

PRINTER IS 1 

BEEP 

PRINT "SELECT TUBE DIA TYPE:" 

1 1 ds “2 

PRINT "1 SMALL" 

PRINT "2 MEDIUM (DEFAULT)" 

PRINT "3 LARGE" 

INPUT I ids 

DDTMTCQ TC f7»1 



in file "" ,10A" ;D_ 



D 1 * . 01 27 ^ 

Hq- ^ 0 1 Q00 I 

D 1 * . 0 1 905 
DZ“ . 0 1 5975 

rTncQ'Tc: 

U/ J • V/ I ^ W I ^ 

IF 1 1 d* * 1 THEN 

n , * (?ig\Qc:oc 

Do“.0l27 
Dr* .01 
END IF 

TC T + TUCK! 

XI •mJ lllk_l« 

aoc 

i_/ w • ^ 

Dr* . 02 

CMH TC 

k_ I XI 



ID o^ m*d*ii 

on 



.1 I I LJ X W I ^ 



on « P 



unnea 1 



, 1 : 



1 » »-> ^ 

X w I ly o I . 



OD of unheeted outlet length 



TO BE MODIFIED WHEN KNOWN 



'TO BE MODIFIED WHEN KNOWN 



1 A QQ 

1 -r ^ w 


Rm-Do 


*',_0G( Do/Di >/(2*Kcu) ! Wall 


resistance 


based 


on outside area 


1 cm 

1 wJ V/ t 


BEEP 












1 C(AA 

{ ^ T 


INPUT 


"ENTER PRESSURE CONDITION 


(0-U,l-A)" 






1510 


BEEP 












1 043 


DD T MT 

1 1 \ X 1 1 


USING "16X, 


, ""This analysis 


includes 


end“f in 


effect 


1 040 


PRINT 


IJC T M0 *' 1 6X 


, " "Thermal condu 


c t i V 1 1 y 


» » •• 


,3D.D,'"' (W/n.K) ;Kcu 


1 54 9 


PRINT 


USING " 1 6X 


, " " Inside diamet 


er , Di 


« •• " 


.DD.DD,"" < nn ) ;Di-'‘1000 
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1552 PRINT USING “ ! EX , “"Out £ ids dianeter, Do 
1E00 IF Ijob-1 THEN 1649 

t C(7lQ QCCD 

I cy W ^ U. ^ J 

1 C 3 C IK*PUT "GIUE ft NftME FOR WftLL TEMPERRiTURE FILE 
IE39 CREftTE BDftT Utf$,S 

1B42 ftSSIGN* @Fils1 TO U*tf$ 

IS49 IF Ijcb-I THEN 
1 C C 1 I C V * 2 
1G54 GOTO 1702 

ICC-? CMH TC 

I w wr I I ^ uy XI 



•' “ ,DO.DD, 



Ut f$ 



|cc0 beep 



1 cct; 

1 «_/ w/ 


T MpiJT 


" select 


OUTPUT ( 1 “SHORT 


2=lONG 3*RftW DftTft 


> 

0 


1 CC Q 


IF If 


TUCM 

^ W 1 1 1 L_ 1 <1 










1 0*72 


PR I NT 


USING 


1 CV 

1 w / \ . 


1 ■ / H 


; instrunenteo 


SMOOTH 


1 C 75 


ELSE 












1 E7S 


PR I NT 


USING ** i 


1 CV 

1 Ur /\ . 


, t UA ur Ui ^ |U w 


: INSTRUMENTED 


FINNED 


1591 


RR I NT 


USING " 1 


1 C V 

1 lU / \ , 


,"“Fin pitch. 


width , snd height 


( nn ) : " *' 


7 ^ nn " 


. c^ c. . 


c u 

,1 * 1 











CO /I Chin TC 

L^w-r 0.1 

702 J = 0 



DD.D0.2X 



. n 






7 nn 



774 


. > w ^0 w wo ^ ■ 


! 




777 


J “ J 1 






'70r7i 

1 W w 


T C T ~.c= 

01 0 t 1 


1 JWCM 




'70^ 

1 w -y 


9EER 






7S4 


OUTPUT 


70Q . "AD 


PP40 AC 40 Upc" 


795 


0 U T P U T 


709 ; " ftS 


Sft. " 


797 


T MQl IT 

0 ■ ^ 1 w 1 


"LIKE TO 


CHECl* NG CONCEN'^RftT I ON 


*7 0 Q 


IF J = 1 


THEN 




*702 


OUTPUT 


7r7iP . "AD 

1 W w> A III. 


AC 40 A^4| MRP" 


795 


0 U T R U T 


70Q . " AO 


Sft " 


799 


END IF 






90 1 


BEEP 






OQ\A 

W V< -T 


INPUT 


"CMTCQ ri__ 


OUM.ETEP REftDING" Fn 


307 


ni ITQ 1 IT 

W W 1 1 W 1 


70 Q ♦ " AD 


ACC0 AC62 UR5" 


9 1 0 


ni ITDI IT 

w w • 1 W 1 


ICAO * "AC 

1 V/ u 1 > I ^ 


Sft " 


9 1 3 


CMTCQ 


709;Etp 




c t c 

w 1 Ly 


OUTPUT 


•7(7lQ . " AC 

1 fc/ ^ 1 I 1 W 


Sft " 


9 1 9 


QCCQ 

w w— r 






07 7 


I N P U T 


" CONNECT 


'JOLT ft GE LINE" Ok 


S25 


C M T C p 


7 0 9 ; E V G 1 




92 9 


BEEF 






93 1 


I N R U T 


"DISCONNECT UOLTftGE LINE" OK 


934 


TC 

0 < ty V w 


K. 1 THEN 


940 


BEEP 






943 


T M p M T 


" T MM A[__ T n 


UOLTftGE - TRY ftGftINi " 


0/1 c 

W-T ur 


GOTO 1 


9 1 9 




34 9 


END IF 






952 ■ 


ni ITQI IT 

W W 1 ' W 1 


7 f?l Q . "AC 

1 «y ^ ^ 1 1 >y 


Sft " 


955 ' 


ENTER 


nOiO . 1 1 + 

t W 1 V V 1 Wl 1 > 



! 000 



2X 



112 



1 QCQ 

1 w w w 


OUTPUT 709;” AS SA” 




1 QC 1 

1 W W 1 


ENTER 709;Bar*ip 




!8G2 


Etp-Etp^l .E+G 




]SG7 


OUTPUT 709; "AR AF40 ALE3 


i IDC •• 

V 1 \ 


1873 


Nn”7+Inwt 




187G 


FOR 1-0 TO Nn 




1 879 


OUTRUT 709; "AS SA" 




1 ooo 

1 u/ w <_ 


IF !>7 THEN 




1 QQC 

1 W 'W/ 


Se-0 




t QQQ 


FOR K-1 TO 20 




1 QQ 1 


ENTER 709 ;E 




1894 


£e-Se+E 




1 QQ'7 

1 W 1 


NEXT K 




1 OrTl^Tl 
1 


Enf ( I )-AB5(Se/20 ) 




1 Q03 


Cl cc 




1 Q(7lC 


ENTER 709 ;E 




1912 


f ( T )mf\OQ ( r ) 




1915 


CMH TC 

%lJ 4. I 




1 Q 1 c 

1 <_l 1 LJ 


Er,f( I )-Er;f< I )*1 .E+E 




1919 


MCVT T 

1 <»l_ /\ 1 * 




1 92 1 


OUTPUT 709; "AS SA" 




1 224 


OUTPUT 713;"T1R2E" 




1 927 


UATT 9 

W 1 I A t 




1 930 


CMTCp 7 1 3 ; T 1 1 




1 233 


OUTPUT 713*"T2R2E" 




1 935 


WAIT 2 




1 939 


ENTER 7 1 3 ; T2 




1 Q 42 


OUTPUT 715;"T1P2E" 




1 245 


WAIT 2 




1 QA Q 

1 w> -r w 


c M T c p 7 ] 3 ; J 1 2 




1 25 1 


T 1 ( T 1 1 u. T 1 2 ) * . 5 




1 254 


OUTPUT 713;"T3R2E" 




1 950 


BEEP 




1 QC 1 

1 ^ • 


INPUT "ENTER PRESSURE GAUGE READING 


1 952 


Dv,o,-> 1 *CQQ/! n 1 pc T rn 

1 1\_/ 


Ps 


1 972 


Cl cc 

l_ k_ ^ k_ 




1 QQyl 

1 v/ W-T 


CMTCD «C-i1«.Dwr>l D4.^ 


= C+^ 

Qil jL.il! 


1 987 


IF J“1 OR J“20 OR J-Nrur. 


THEN 


1 QQITI 

1 ^ XJ 


Klr^ae 1 

•'•a ' 




1 99”^ 


ELSE 




1 QQC 

1 w w 


Mr>,sr7i 




1 999 


END IF 




2002 


CMH TC 

^ 1 ^ XI 




2009 


T 5 1 csn 1 ■FNT V 5 v57( Er; f ( 0 ) ) 




r>CA(AQ 

t- x/x^ ^ 


T5tsar}l“Tstc5n1”“273. 15 




201 1 


T 5 1 esr,2 *FNT V5 v5S ( Er. f ( 1 ) ) 




9^71 1 9 

<- V/ 1 i- 


Ts tsar.2'T5tean2“273 . 15 




\ A \ 

4- »./ t -r ; 


Tstcan-( Tatcaml +T5tear;2 )/2 . 


2015 


Tstaar.-Tstsaml 




9^71 1 9 

4 . • 1 


TrooPi-FNTvav58( Enf (2 ) ) 





{ Dr>. 

' ' a'3 



P a 

,r 






vsp ] 
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20 ] S 


Tccn“FNTvs v59 ( Enf ( 7 ) ) 


20 1 9 


Trccm“Troon~273 . 15 


9^71 '>r?i 

4- ^ C, ^ 


t^qj-j«t^qj-;_273 , 13 


202 1 


t uu r 1 xj 9 


2022 


FOR >0 TO Inwt -1 


2023 


Tu;( I )*F.NTvsv57( Enf ( 1-^-9 ) ) 


A 

-r 


Tt4( I )*Tu;( I )-273 . 1 5 


202 S 


Tu;r.=Tuin'i-Tw( I ) 


7r7i'7 O 

4- ^ 4- vJ 


NEXT I 


2032 


Tojr^“Tu;!^/Inu-t 


2044 


pe*4»pMpw«+( ) 


2057 


nilTPHT 7r7iQ."AD ACC/! Al C /} 1 JDC " 

wwtl WJ IM\ III U/^ iJU-U/-^ 


2059 


OUTPUT 709* "RS Si9” ^ PRESSURE TRANSDUCER 


r>(hzi □ 

^ tu ^ -u 


S s“0 


r>(hC,CA 

4— XJ \J 


pQp TQ 20 


2092 


pMTpp 709;£trfin 


9/TlCC 

4_ W 


c«=52-fptr3n 


2099 


NEXT K' 


2097 


Ptr3n*^ABS( Ss -^20 ) 


207 1 


BEEP 


2072 * 


PRESSUPE IN Ps FROM TRANSDUCER 


20’^3 


Pvs'^2“( - 2 . 93904*Ptr3n-^l 4 .7927 '**9994 . 7 


2079 


D+ 4 a« 4 . 1 =Dw=»r^ 1 1 CAllCir 

iVwWbi WtIWWi— 


2077 


Pts5t2“Pvap2 ^ TRANSDUCER 


2090 ‘ 


CORRECTION FOP AIR CONTAMINATION 


209 1 


DL P=Py op 2 * ] . 3-3 1 TRANSDUCER IN hPa 


20 CT 


pL «=p==j+*i ,p -3 1 SAT. ^PESSU^E IN I'Pe 


"•'034 


FF g’Pvap 1 * 1 . E~3 ' GAUGE IN !■ 


2092 


T 2 a + = p M T 2 p ( P S 5 t ' 


2093 


T 1 =CMTv.c m ; Dv, =u,A 1 

V p , . 1 X 1 . p , y ws p 1 


2034 


T V . = 9 = C M T V . c rx ;• D w =» n ' 

IV(_.P4_ i.yivop-. 


903 q 


u 5 1 = F N U 5 1 { T s t e a r, ) 


2 1 I 0 


N u; “ 1 9 015 


9 11 a 

2119 


IF I f t ” 1 THEN M u: V *= 1 9 7 . 3 9 
IF Ift = 2 THEN M’uu'v=92.05 


21:7 


M^.r,r^l3ir7ir7i*/D.,«^1-Deaf\//Dv,afAl-/t _/M..vy/9Q aC^^♦D«a4•^ 

\IYLJipi lvJOt4/V|VLjpi ^l» \lllA V/l_W.^l-»f4 IwJdv/ 


2119 


1 / Dw»n'^-C3e ^ Dw = rx9-/ 1 _/M,,w/9Q qc\N^Des»+N 

11111^4. . XJ XJ ^IVWp4_ l_/^<4//Vlvdp4_ \i< \IiUtv/^w.duj// iddc/ 


2 1 34 


IP Icv“2 THEN 


2135 

9 1 ■77 


PRINT 
P R I N T 


9140 


DDTMT MCTMC " 1 Y " " R a 4 a «= «» 4 ” HR " . T 

ll'lixi WOXIXW 1 XJ t\ f U/U1«>W U'wb IM.^IlU'wl }U/U/ 1U 


9 1 4 a 


OUTPUT 709; ‘'AR AF40 AL40 UR5" 


2149 


OUTPUT 709; "AS SA" 


2149 


CMRTC 

U- 1 X U/ X 1 • 


2150 


PF I NT 


2 1 52 1 


TC T^w«9 AMR TUCM 

XI xwy 4_ IIIHL^ ' Ilit^lX 


2 1 cc 


PRINT USING "IX "" Psat Pga Ptran Ttrans Tsat(P) NG1% NG2% 


2 1 59 


ppT^vjT ‘'SING "IX "" ( kPa ) (KPa) ( kPa ) (C) (C) Mcial Mo la 
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21B1 PRINT USING " ! X ,3( 30 . DO ,2X > ,2 ( 3D . DD ,2X ) ,2 ( M30 . D ,2X ) " ; Pks ,Pkg ,Pkp ,Tvap2 ,Tsa 
t , M f ng 1 , M f ng2 



21S4 

21G7I 

9 1 noi 

<_ I t ^ 

2173 

'P 1 7C 

I I W 

2 1 79 
2 ) 82 
2 1 SB 
2199 
2191 

2 1 94 
2 1 97 

^ Wr W 
4. ^ ^ W 
4_ 4_ V/ LJ 

2209 
22 1 2 

22 1 5 

22 1 9 

222 1 
2224 

r> 7 *7 
'P9“3:<7i 

i_ <- w ^ 

2233 

2236 

7 7 A Q 

2242 ■ 

2245 

2245 

2247 

2249 

2249 

77C<7l 

^ ^ w 

2251 

2252 

2254 

2255 
2255 
2257 

9 7CQ 



THEN 
AND Ng-1 THEN 



PRINT 
END IF 
IF Mfngl>.5 

DCTD 

w ^ I— > 

IF In-1 
BEEP 
PRINT 

PRINT USING " 1 0X , “ "Encrgi zc the vecLiun eysten 

BEEP 

INPUT “OK TG ACCEPT THIS RUN ( 1 -Y ,0=N )? '* ,0k 
IF Ok-0 THEN 

DCCD 

1^ I. I 

DISP "NOTE: THIS DATA SET WILL BE DISCARDEDi! " 

WAIT B 
GOTO !7S0 
END IF 
END IF 

TKin TC 

I. I ^ XI 

TC T«»1 TUCM 

TC c*M./|0 Qp p~',>]00 THEN 

I f r,-0 

BEEP 

TMDIIT " TMrnDDCPT CM M«A7rCDT « HC I C T C R C C AM I T ^ ^ ” T ^ 

XIO) Wl XI'IWWI\l\k_wl I II \ I nwwt_l t jW' L^b.l_t_lU.%Ld<l_l ilWl—l / / ,1)11 

TC T^)V»«:<7I TUCK! 1QR/I 

XI XIII 1-' llll.lll IWV/*^ 

CMR TC 

l_ I H I-/ XI 

CMR TC 

l_ I XI 

AKiA|_YSIS BEGINS 

T , 1 5= C M J w e w C O ( p jv; f ( 3 ) ) 

T 1 9„CMTx^.e wCC ( Cj^ f / C ) > 

FNTv; 

CMTw ^ ^ ^ 

T , 1-997 



1 

Tcl 

Tc2 

T 1 t 



> V 5 9 ( E n f ( 4 ) ) 

: v.CC/ { c \ \ 



j w \ iz. n I 
1 c 



.T ,9-997 



! C 



T 1 2 " 

Tcl* 

Tc2*Tc2' 
Tdei 1 -Tc I 
Tdel2=Tc2- 



1 5 
15 



f(3 )+Etp/20. 



27 

273.15 

1 -T i ] 

Ti 2 

Tde 1 3=T2“T 1 

Cl 1 =c 

I- K H ' •- 

n + . 593 1 E-2 - 1 

A 

2259 Tne-Dtdc*Etp/1 0 

2262 To3*Til+Tne 

2263 PRINT USING 

RCI T9: TD TI C " » 

(xr l_ l_ I W II X l_ t_ 

2264 PRINT USING 
2255 PRINT USING 



C- 1 Q*C+ 1 

l_ I W i_ ». p I 



1 V 



1 V 



50464£-5^Etp] + l . 2 1 70 1 E- I 0*Etp r 2-5 . 1 1 54E- 1 5>Et p r'3-f 3 . 220 1 



TINl TOUTl TIN2 T0UT2 TIN3 T0UT3 DELTl DELT2 



( Teflon ) 



(Meta] 



( Quartz ) 



99CC 



C<- 1 mADC f T 1 



"2X,10(2D.DD,2X)“;Til ,Tol ,Ti2,Tc2,Tl ,T2,Tdeil ,Tdel2 ,TdeI3 ,Tn e 
TI ) 



115 



PPC'? 

4- <- W • 


Er3"A5S ( T 1 2-T 1 ) 






2259 


PRINTER IS 1 






2270 


T p Er 1 > . 5 OR Er3> 


.5 AND In=] 


THEN 


227 1 


BEEP 






2272 


PRINT "*^CT AND TC 


RTppcp pY 


MRDc THAN 0.5 C" 


2275 


Ok 1 = 1 






2279 


cTMn Tcr 

U. 1 t XI 






2297 


jp R[.]*0 AMQ Erl> 


.5 AND Ir*1 


THEN 1790 


2299 


IF Ok 1-0 AND Er3> 


. 5 AND I n“ 1 


TUCM ] 700 


2290 


Er2"ABS ( ( T2-T 1 )-( 


Tris ) )/(T2- 


T1 ) 


2293 


TCr C^'^\ RC AMR T «A 

Xi L-I i_.- .ow XI'I 


=1 THEN 




2295 


BEEP 






2 2 99 


print "QCT AND T- 


PILE DIFFER 


BY MORE THAN 5%" 


2302 


0 K 2 * 1 






2305 


TT RLT«R AMR C^r>\ 

XI wi-4_ XU iiiykw^ k_i X..' 


.05 AND In= 


1 THEN 1790 


2309 


END IF 






2311 


pPJMTpp TC 70] 






2 312 

1 1 ~'i 


INPUT "COOLANT TEMP. RISE ME AS. (0“TEFLON 1*=METAL 


2314 


Tp JtR*0 THEN 






2 315 


C D T MT 1 IC T MR " 1 V " 

I I > X Ml I W W X 1 . /> , 


"USING SINGLE TEFLON THERMOCOUPLE 



2317 

2320 

9 "?, *? "7, 

2326 

2327 

o ^ p 

'^3'^ S 

2335 
233S 
2335 
234 1 
'’344 

2347 

2"^ 4^' 

'? T. A a 

2350 

235 1 

2352 

2353 

2354 

2355 
2359 
2352 



T 1 i = T , 1 

T 2 c “ T c 1 

CMR TCr 

L_ 1 > W XI 

T p T + I TUCM 

PR I NT US I NG " 1 X 

T cc 

T 1 i 1 2 

T2o* To2 
END IF 

ip T + jyv = 2 THEN 
PRINT USING "IX 
T 1 1 “T 1 
T2 c*T2 
END IF 

IP Itr^"3 THEN 
PRINT USING "IX 



T 1 
T9 



T 1 ! 

T c 3 

CKin TC 

l_ I <1 W XI 

T ^ T 1 , 

I w y y > 11 , 

I f t * 0 

Pr^l 



,/ T . 



' y " 



•IICTMC CTMRI C MCTAI CUtrATUCm T U CT D M H P R I I C CRD RRRI AMT TCM 

\_/,,/xiy^ wXiywL,.!,. iit_i ii_ ^iii_iiiiii_u/ iiii_i\iiOwwwi i_i_ i wi> \^W>./L_Miyi ii_ii 



'USING ^UHRTZ THERNONETER FOR COOLANT TEMP. RISE""' 



'USING THERMOPILE FOR COOLANT TEMP. RISE 






iwuc I i«ii\Iiwuu \ luav^ 

Md*(G.7405*Frv^l3.027 )/1000. 

Md=Md*( 1 .0355-1 .95544E-3^T1 i + 5 . Z52E-B^T 1 i"2 )/l .0037 
Mf-Md/Rhou) 

VIXI li|/>IX U> J. ^ t ~T t 



2355 
2359 

2371 

2393 IF Inn-0 AND 'Jw>.5 THEN T2c=T2o- ( -2 . 73E-4+1 . 75E-4^'Jw + 9. 35E-4^|'v'w''2- 1 . 96E-E 
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U. . ^ A \ 

V WW w/ / 






2389 


IF Inn-! 


THEN 


2407 


IF Inn-2 


THEN 


2409 


IF Inn-3 


THEN 


2413 


Q-Md*Cpui<‘ 


( T2o 


2415 


Qp-0/(PI* 


Dc»L 



o-T2o“( -E . 44E-5+ 1.71 E“3 ^‘Juj+4 . 45E“4^Ulj" 2 + 4 . 07E“5^‘Ju' 
IF Inn-2 THEN T2o-T2c-( -3 . 99E-4 + 2 . 75E-3*Uw+ 1 . 4SE-3*Uw"2+9 . 1 6E-5^U-^' 
T2o-T2c- 



24171 

2419 

2419 

O/t 

4- -T 4_ W 

242 1 
2422 



ITERATE TO 

Tj.f M 1 ■» T 



FIND OUTER WALL 

un } - 1 wn 
Twr.2-Twn 

Turn 1 - ( Tur. 1 +Tu)d2 )/2 . 

Ct-0^L0G( Do/Dr )/ ( 2 . ^PI ^Kcu^L ) 
Tur}2-Tu;n+Ct 



TEMPERATURES 



0 -<1 9 7. 
i_ -r 


IF ABS(Tu;n2-Tun1 


)>.001 THEN GOTO 2420 


2424 * 


ITERATE TO FIND 


INNER WALL TEMPERATURES 


2425 


X 1 I ^ ^ ic X I ,f ^ 
1 wr 1 ^ 1 1 < 




2425 


Tun4 = Tu;n 




2427 


Tu;r;3 = ( Twn3+Tujn4 ) 


/2 . 


2429 


P4- , =n*i nt7(' n.r. /n , \ 

w L~wwvLri / uyx/ 


/ ( 2 . p I * K c u * L ) 


242 9 


Tu;n4=Tu;n“Ct 1 




9 /I 

i. -r ^ 


If ADC ( T^^fy»4_Tj^p«3 


)>.001 THEN GOTO 2427 


2431 


T uj r 0 — 0 




n /I •7'7 

A- “T ^ 4_ 


T. , a/7l 

< ww I I X w • 




9 A 77 

4- -T ^ ^ 


CHD T*n TO T^,.,+ - 

1 sJl' ± XJ IW 


1 


2434 


T. T \«T. .( T ^x^'4 

IIAiVj. /> XJ 4, 




9 >1 7C 

4_ -r _i 


T. . W T T )-r+ , 

IIA/XVJ./ / t X 




947 c 

1_ "T 0 


1 WI'IW 1 U/l IW ■ 1 lA/W \ X / 




2437 


TLjri-Tu;r)i-^Tu;i( I ) 




9 A 7Q 

-T W 


N EXT I 




9 4 "T** Q 


T jv, Q ss T fvj ^ / T + 




2440 


Tun i '==Tun : / I nu t 




9441 


PRINT 




2 443 


PRINT USING “IX 


1 w^xvXw}i 




c 

XJ 




9 /j 4/1 


DD T MT MC T MrC *' 1 V 

' ' ' X .«IW . /V , 


^ 9 \ i.ial 1 

wL^Wwi 



f ncc n \ 



•• •' ,B(DD.DD, 



1 y V : 



) 

2445 
2447 
c Tu;n 
2449 

9 /I Cf7i 

4- “T wj 

245 1 

2452 

2453 

2454 
2457 
2459 
2459 

ACa 

i_ -r w w 

240 1 
2452 I 



PR I NT 

PRINT USIN6 
PRINT 

T c + «r7i 



OLilw^ 



snd inner wai! 



A rfs Q 



,2(DD.DD,1X )" iTuir 



•FNK 



. T awr. \ 



Muu-FNMuu;( Tavg > 

M, .... . =CMM... W T,.,,v, , \ 

I I'^IIULVJV IWfiX / 

Rai-Rhcw^^Uu-^Di/Muu) 

Prw-FNPru!(Tavg ) 

Dtu-Q^^L06(Do/D: )/( 2*PI^Kcu*L )^.5 
Lntd-(T2o-T) i )/L0G( ( Tuni-T ! i )/( Twr;i-T2o ) ) 
Perin-PI ^‘D: 



C 1 1 ^ ^ 

w u t I Q j 



“P I *D i 



PRINT Rci ,Prw ,Dtu; . 



<4 A 0^.~ i , 



curfai .Q 
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t>J OJ 



I 



24B7 


ftnarea!-PI‘(Dc''2-0r2 )/4. 






24B8 


ftnarea2-PI»( Dr''2-Di''2 >/4 . 






2 4B B 


H,i- 0 /(SLirf 3 i*Lntd) 






2470 


Ff l-(Kcu^‘Pcrin*ftnarea! )'• .5 






247 1 


Ff 2 ”(Kcu*Pcrir',*ftn 3 rea 2 )'' .5 






2472 ^ 


PRINT ftnarca! ,ftr.area2,Hi ,Ff1 ,Ff2 






2474 


u . t _u , 

f 1 X 1 1 1 X 






2475 


H 1 Z H 1 1 






2475 


{Hi 2 -Psrin)/(Kci;*ftn 3 rea 1 ) )" 


_ c 




2477 


P.i3r.2-( {Hi 2 *Perin )/(Kcu*Ar.5rc32 ) )" 


.5 




2479 


11^ lit ^ \ i tAtl 1 A til/ 






2479 


Pf4 = Pf2*Hi2'' .5*(Tumi-T2c ) 






2490 


pLinc5=Pf 3^PNTanh^ AlarjMLl ) 






'? A C 1 


Ftjr.c 7 *Ff 4 *FNT 3 nh(Aia!^ 2 *L 2 ) 






2452 


ir.,^r-.0«U, . * { Ti.ifv,, -Tav,r^ \ 

1 l^liww 11X4. w/l.lt lUiX <- IlMlix 






2493 


FLinc,^. = FDncB-^Func7 + Func9-Q 






2494 


D f Ljinc 5 * ( . 5 ^ F u nc 5 / H i 2 ) -^ 2 . * F f 3 * L 1 / ( 1 


XpP| 0 Qc|-j( 2 




2495 


D^unc7®( .5*FLiric7/H i2)'^2.*Ff4'^L2/^v 1 


xCMPr, « k / P 
• '1 




2495 


DfLinc2 = SL;rf 3 i*( Tuni-Tavg ) 






2497 


O^ur.c/ =DfuncB + DfLjnc7-^Df unc9 






2499 


Hi 1 =Hi2”( FLAp.C " /DfLir!cx ) 






2499 


Tp Pigc ( H 1 1 -H i 2 ) . 05 THEN SOTO 2475 






2 4 90 


H 1 = H i 1 






C i7i 1 

U' U- 


DD T KiT U , 

1 1 \ X I t > IX 






252 ^ 


iliy 






2527 


Tfi lp*T5tcsp/3 + Tu;pc'^2.^3 






2C"70 


p f ssCMp ^ ( T f ^ 1 jvv ^ 






2-33 


[j ^ = P M D ^ 1 • I T » 1 tv \ 






25”4 


Tr- 4 : = U MPr^i . Z- T P , 1 tv . 

WJ_ ilSVrfplA/\l'.X.i 






”'735 


= TP, 






^c~?2 


r. j =P00 






2 S'^B 


riiQcc=ro/( 100-2C )*(T5+»ajv%~Tf-3Q,vj ; 


1 TO DC MOnTCTCn 

<W WL. XU»W 


2575 


rp„^s-pMPp,^( ^ Tccp"^'''5t33n .5 ) 






2534 


jVl r4 w = 0 






2597 


9 p = ( 5 '/ 0 1 * ^00 ) '' 2 / 5 . 7 5 






*^590 


U«r-=rr-.c^x 'i 

iiw/^ l.tt.W.4,, l>^OIl. 







2593 


lli^vw V \ X^p UtX^w/.t / ltl.tv 


U« ^ N /UPr^ 

1 1 vt w / r III y 


2595 


T p Pi g 9 ( ( N d V N d V c ) / M d V 


c ) 7 . 0 1 THEN 


2599 






2502 


GOTO 2593 




2505 


CMR jp 




2509 


N d V = ( N d - N d c ) * .5 




'^C 1 1 


Ut^ = CMIkt«i^TeP-..«, \ 

vy l'«V»./VV l«JWWUlll^ 




25 1 4 


Uvt = Mkv,*llr^ / Av 

V » . I w V V y / 1 , .1 




2524 


Tr4r^P«T«» + -T. .fwt-K 
■ i.t w 1 1 w 1 vui'iw 




2534 


Hc*Qp/Tdc f 




2535 


( /i/ ( 

ii 




2535 


R c t p * ( U V * D c ♦ Px h c f ) / M u f 




'T>C-'7 

4 . 1.1 ^ 1 


M.j . « » { * Hq ) / p f 




2539 


Nurcl-Nuss/Rstp" .5 
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2£39 Ff-9.91 *Dc*Muf^Hfg/( Vv"2^Kf ^Tdcf ) 

2640 Hnu55-.728-(Kf "3^9.91^Hfg^Rhor2/(Muf*Do^Tdcf ) )".25 

2642 PRINT 

2543 PRINT USING “ 1 X , " Uinf Tdcf Ofl Ho 

2644 PRINT USING “ 1 X ,2( DD . DD ,2X ) ,2( N2 . 30E ,2X ) " ; Uv ,Tdc f ,Qp ,Hc 

2645 PRINT 

2546 RRINT USING “ 1 X Nj Rctp NuRc Pr 



inu 



2649 PRINT USING “ 1 X ,2( NZ . 3DE ,2X ) ,2(DD.DD,2X ) ,2(MZ.3DE,2X )" ;Nlis 3 ,Rctp ,Nurel ,Pr 
Ff ,Hrrjss 



9CC/1 

1 - w taJ -r 


Rat -Now/Nuwi 




2064 


Ar^. .e «u ^ ^ /u. . 

iia’^L^a / i\w 




2674 


Ar.» 1 l. .«n 1 A.P kr-^.., 

III \\/WU U>A IXItWIAl 


)/Nou] 


9C'7C 

A. W 1 


Piprs*Nou;"^Cpu;/Ku; 




9C9C 

^ W 1 w 


Stco-iPnus/CPiras" 


. 9*Pipr 


2677 


Nupr^Piros / ( Aprs 


. 3333^ 



3''.3333-^R3t^’.14) 
'vct" . 14 ) 



2679 PRINT 

9CQ/7i DDTMT riCTMC "IV 



Ucoci 



Rec 



i tcc 



200 1 DDTMT mctmc "]v 0(M7^3nc 2X)”;Uui Rsi Hi ^^.nus Stcc Nupr 

2692 INPUT "CH.^iNGE TCOOL RISE'^ 1 2“N" lir 

2693 IF Itr“1 THEN GOTO 2312 

2694 OUTPUT @F i i e 1 ; T-oJC ( * ) 

2704 IF Ir;=1 THEN 

2707 BEEP 

2710 INPUT “OK TO STORE THIS OATi^ SET ( 1 “Y 0“N)'?'' Oks 

7'71*7. TC nL»«1 TUDM 



2725 OUTPUT 0File;Bvcl Bar.p Ptrc 

2739 0MH If 
DDDD 



D+r. c:.v^ T1 T7 Dv.r.r.1 



2761 INPUT “UIILL THERE BE ANOTHER RUN ( 1 =Y O^N)"?" Igo 
2704 Nrur,*J 

2767 IF Igo=1 THEN 1777 
2770 ELSE 

2773 If THEN 1777 

2793 ENO IF 
2933 BEEP 

2S36 ‘ INPUT “ENTER PLOT FILE NPiNE" FplctS 
2Q30i -ASSIGN @FiIe4 TO F"^Iot$ 

2942! FOR I“1 TO Nrur. 

29401 enter 0Fi 1c4;G'^ Ho 

29001 M0VT I 

^ZC7\fh^ DDTKJT IICTMC '* 1 fTlV " " MOTCT • "" 77 

n ^ ^ 1 

i-r » A A C 

3007 BEEP 

TC Tn^-1 TUCM 

^wix^! Ai X|J 



data r'jDs uere stored in file "“,10A“;J 



3013 PRINT 

3010! ppjMT ncTMg *' 1 0X ““NOTE* 
“ “ ,101^" ;J ,P_fiie$ 

■7.(7\^ Q ] CMH TC 

I ^ ’ 1.1 <tL^ A I 



were stored in plot data fils 
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3073 


ACC TCM «c: ; 1 ^ TO ^ 

\rt xaC «w ^ 


■7f7k'7C- 

^ W 1 U( 


ASSIGN @Fiic) TO ^ 


3079 ! 


ASSIGN 0Fiie4 TO - 


3092 


IF Iso-2 THEN CALL 


3092 


IF Iso-3 THEN CALL 


3)15 


END 


3)19 


ncp F^'Pvst(Tc> 


3)21 


CON /Fid/ Ift 


3)24 


DIN K ( 9 ) 


3 1 27 


If Ift-0 THEN 


3)30 


HAjA - 7 . £91234554 


3 ) 33 


HATA 4 .) 57 1)732 20 


3)35 


READ K ( ) 


3)39 


t»/t^j.273. )£ )/E47. 


3142 


C « 0 


3)45 


PQP ^=0 TO 4 


3 ) 49 


S l; n - S 0 r + K ( N ) ^ ( ) - T ) 


3)5) 


Kj 5 V T M 


3)54 


D - = C .J fv^ / / T ♦ ( 1 X p ( 5 ) ^ 


3 1 57 


Pr=EXP ( 5r ) 


3)50 


0 = 9 0 1 r)c?taaa*.D^ 

1 <_ 4_ 1 4_ W W 1 1 


■7 1 CT 

1 u w 


END IF 


3)55 


ip Tf+c*i then 


3)59 


t^-Tq*.) ,94-32 + 459.5 


3 ) 72 


D = 1 a'' / y.‘7. A-^.y.a 

1 t XJ ' w/ . W W ^ ^ 'U 


3 ) 75 


p»p* 101 329/ 1 4 . G9£ 


3 ) 79 


END IF 


T 1 O 1 


IF Then 


3)94 


A - 9 . 3 9 4 5 9 5 - 3 0 5 5 . ) / 


3. 1 57 


P=)33.32*)0'’A 


3 ) 90 


END IF 



3 I 93 
3 1 9E 
3199 
3202 
3205 
3209 
32 1 ] 
"^214 
32 1 7 
3220 
3223 

w 

3229 
32'^2 
3235 
3239 
'^'^4 1 
3244 

/I '7 



T \ / M4. 1 ^ 



-7T 1 C \ 



DCTI IDM 

i > L_ I W I ^ I '< 

CMCMH 



nrr t \ 



QQM /Pi|j/ + 

T p I f t “0 THEN 

= / T-1/71 S 

lll^ \t iiy/ 

END IF 

IP I f t * 1 THEN 

T f »T* 1 ^ Q + 32 

Hfg-7.0557957E+1 “Tf *( 4.S39052E-2+1 . 25 1 904SE-4^T f ) 
Hfg-Hfg*232G 

cr M n T c 
i_ I X 1 .^ j. r 

ir Tf + o:7 TUCM 

J.ib Ifiu-ix 

T t e.T-f-273 , 1C 
u;^rx*i / cxc_T t . 

Ill^ I •W./^U'Tk— *U/ II' 

END IF 
RETURN Hfg 

CMC Kin 

I I 4 ^ I X i_i 

ncp pMMuw(T) 



120 



COM /Fid/ Ift 

TC T f ♦ _pi TUCM 

XI XIV \J IIIU.IH 

a -247 . 8 / ( T+ 1 33 . 15 ) 

Mu-2.4E-5*I0"A 
END IF 

IF Ift-1 THEN 
Mu*9 . 962 98 1 9E“4“T* ( 1 
END IF 

IF Ift-2 THEN 
Tk-l/(T+273. 15 ) 

Nu-EXP(-l 1 .0179+Tk*( 1 . 744E+3-TK ^ ( 2 . 90335E+5-Tk ^ 1 
END IF 

DCTIIDM M.. 

)Vl_ i Wl\l<l !»«_« 



1094608E“5-T*5.566829E“9 ) 



! 265 1 E +9 ) ) ) 



ncc CMllx,e + ( T 4 - 

I-/ I— I « i<<s/vov \ I V 

QQM /P i J f t 

IF lft*0 THEN 

DcCMDv/ e 4 / T 4 . N 



’ 1 . E-4 ) 

X ( 1 - c y p ( _ V ) ) 2 . 5 

r7K7MC*Cl_ f7»(7H7tQ/l7 





1 Cl' 7 ' 79 C 9 _T^*C 



3355 


CMR TC 

^ 1 <1 U/ XI 


3359 


RETURN U 


aTC 1 

W w/ W 1 


FNEND 


3354 


ncC CKlPr.. T N 

i^U.1 1 i<iwpw \ 1 / 


’7TC'7 

^ W 1 


roM /r 1 r 4 / T P 4. 

W^ll / IXW/ XI V 


A *7 (7k 


TC TP 4-= (7k TUCM 

XI XIV lllk_l<l 


^ ^ 1 w-/ 


r^s. ,^A 7 11 7r7kQCQ 

wpvw *T.4_I I4.VWWW 


3375 


END IF 


3379 


IF Ift-1 THEN 


3392 


Cpw“9.2507273E 


3395 


END IF 


^' 7.00 

w ^ w w 


IF Ift-2 THEN 


339 1 


TK-T-^273. 1 5 


3384 


Cpui=4. 1959^( 1 . 


-X-Q -7 

^ w/ 1 


END IF 



.423BlE-5'i'2.71 429E”7*T ' 
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3400 

"a ^ ”7 
^ ^ 

* 7 y 1 /Tie 
^ -r ^ O 

3409 
34 1 2 
3415 

A>1 1 Q 

W/-T I W 

3421 

-7424 

3427 

3430 

3433 

343c 

3439 

3442 

3445 

'T.A AQ 
W< -r -r w 

345 1 

3454 
3457 
3450 
3453 

3455 

345 9 
3472 
3475 
3479 
3491 

3404 
“ u'7 



DCTllDM P.A.,* 



Ul* 1 000 

CMCMH 

) I I) I 0 w 

DEF FNRhcw(T) 

POM /U ^ A / T ^ + 

W\_/ll / < At_»/ 111. 

If Tf+*0 TUC|\j 

Rc“999.52945 + T^( . 0 1 269-T* ( 5 . 4925 1 3E-3-T * 1 .234147E-5 ) ) 

END IF 

IF Ift»l THEN 

Ro*1 .5207479E-^3-T-^( 2.21 96346-^1*2. 3579291E-3 ) 

END IF 

Tc ifi^2 THEN 
ti.bTx233^ 15-339. 15 

Mf«9 24949E-4-?-TK*(5.2795E-7^'Tk^( 9.2444E~10+Tk*3.057E“12 ) > 

p Q B I / n f 

END IF 
RETURN Re 

p M p M n 

npp p}vjpp.^(T) 

p p - p M P p .J ( T ) * P Kj M ^ ( T ) / P M 1/ ( T ) 

RETURN Pruj 

PMCMn 

t I «i_ I « u/ 

n r p p M K ( y ) 

PPM /f\^/ I f t 

IF I f t *0 THEN 
y»(t^273. 19)/ 273. 15 

L',. = _ Q97 /I 'Tj- Y * ^ 7 p7:QC:-y*/l OPPi 7_Y*/ C9C'7'7- Pl'7-5:4./l * V ^ n \ 

PNin TC 

I . j. r 

IF I f t * 1 THEN 
Kuj-9 .2095239E-Z-T^ 

CKin ir 



* ] 429PP~4-^y * 2 . 3909524E~9 ) 



■7 /I QP 

w. -r ,_J w 


IF I f t **2 THE! 


3493 


tl=jx273. 19 


3495 


p^= 4 ^ iocQr _4 


3499 


END IF 


“5:cp7 

i- 


return Ku; 


3505 


FN^ND 


•7CPiC 

wJ ^ W W 


DEF FNTaph(X 


35 1 1 


p«rvp( Y ) 


-c 1 4 


n = CYp( ^Y ) 


35 1 7 


T D-n W M 

iwitit SJt ! ! >1 


3520 


RETURN Tape 


3523 


F N E N D 


3525 


npp pMTwoy(M 


352 9 


CON /Cc / C ( 7 


3532 


T“C ( 0 ) 


3535 


FOR 1=1 TO 7 


3539 


T«yxP( T )*u*^ J 


354 1 


NEXT I 


3544 ! 


T-T-H4.73385E 


3547 


RETURN T 
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"7 ^r?! ”7 

w I ^ 

37 1 3 
3723 

"7 "7 "7 "7 



FNEND 

QCp pMU f ( T ) 

con /Fid/ Ift 
IF Ift-0 THEN 

Hf-T*< 4. 20394 9-T*( 5.891 32E-4-T* 4. 55 16031 7E-6 ) ) 
END IF 

IF Ift“l THEN 
T f *T^ 1 . 8+52 

uf-Q^2079571+Tf^( . 19467957+Tf*1 .3214296E-4) 
u f -p{f ^2 . 325 
END IF 

IF Ift-2 THEN 

Hf*250 ! TO BE VERIFIED 

END IF 

DCTTI IDM U f * 1 aP\Oi 

ill 9 XJ 

FNEND 

nrc rMTw- / p \ 



C 1 *7 

W 1 


T. .« 1 Q^ 7 > 

1 U> 1 ^ Vj* 


35 1 £ 


T 1 * 1 0 


- 7 ,C 1 □ 

^ LJ t wi 


Ts“( Tu + Tl )* . 5 


W LJ i. 1 . 


.CKiDx,* / T« \ 

tl^ivoyviuj/ 


3E25 


T p 4 ;DC ( ( P_Pq ) /p 


■ 7 C 7 Q 

LJ <_ W 


Tp pq<p then T1 


363 1 


jp Pc>P then Tl: 


3534 


GOTO 35 ! 9 


3537 


END IF 


3540 


RETURN Ts 


- 7 C /I 7 

^ LJ *T w/ 


FNEND 


3E53 


DEF FNTvsv55( U ) 


3553 


con ^Cc55/ T55( 


3573 


T=TCC( 0 \ 


3583 


p Q P Tel T n c 


■!^C Q-!: 


T* T xTCC / T N * 1 1 ^ T 



t ' m I 

= T a 



THEN 



NEXT I 
RETURN T 
FNEND 

nrn: tmt 



' 3 V WLJ XV/ 

rOM /Pr^CC/ TCCrC\ 

WW»I /WW<-^CJ/ 



• 773 c 


TeTCC f r?i ^ 

■ 1 U V / 


T’TaC 

w/ < W«_> 


FOR 1=1 TO 5 


3737 


TeTxTCCi' T N*tJ''T 

1 1 • 1 .-M_» ‘ X / V X 


3747 


M E X T I 


3749 


RETURN T 


374 9 


rKiTKin 

1 < y i_ 1 X !_/ 


3753 


n p p pMT w = w c 7 ( n ) 


3754 


POM /P^C7/ TC7^C\ 

WWII / WWW</ lw>IVW/ 


3755 


T«TC7/ a \ 

I 1 ^ 1 \ w / 


3755 


cno T«1 TP c 

1 W 1 \ JL 1 1 W W 


3757 


T-TXTC7( T ) T 


3757 


MCVT T 

) y k. 1 X 
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37G9 RETURN T 
37 c Q FNEND 

3773 DEF FNTvsvE9(U) 

3774 COM /Cc59/ T59(5) 

•X'lnC. T*TCQ^ Qs \ 

w> I I ^ > I ^ w \ / 

377E FOR 1*1 TO 5 

-xnnn t*t xtcq / t \ t 

3797 NEXT I 
3707 RETURN T 
3007 FNEND 
3917 DEF FNCcshCX) 

3927 P“EXP( X ) 

3937 0*EXP(“X) 

3947 Cosh'.E-'-^P+Q) 

3007 p0Tijpf.j Cash 

3007 FNEND 



3977 
3979 
397 9 
3990 
399 1 

3902 



3984 

3995 

TQQC 

w/ U/ W W 

3907 

3999 
3 959 
3990 

300 1 

"^992 

3993 

T Q Q ^ 

3995 

3900 
3007 
30 QQ 

3999 

3900 

3002 
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Cl ID D . 

w/ W ^ I ' W W 

rnM /c ] M / j ^ 

WWII / I i.W/ 0.1 o 

DIM Emf ( 13 ) 

BEEP 

INPUT "ENTER FILE NFiME" FilsS 

ACC T CM @File TO FiieS 

BEEP 

INPUT "ENTER PRESSURE CONDITION ( 0 =U 
INPUT "ENTER NUMBER OF RUNS" Nrun 
PRINTER IS 701 
PRINT 
PR I NT 

pp T MT ijc T Mc '• 10 V ''"File Ner e • 

TD Tr,^c.r 7 l TUDM 

01 tt.i-.'. 

p p T M T ij 0 T M 0 "10V “ " p jr- « c c ^ « C c r. d 1 1 1 0 n 

Dl CD 



1 *A ^ " T , 



PRINT USING " 1 0 X ""Prt 
END IF 

DKITDD Tr^r^ 

\ri oow^Aiy.oiiii 



iiLjowowji* 



PRINT 

PRINT USING " 10 X,""Dats Uu Tin Tout Ts 

FOR >1 TO Nrun 

ENTER 0 Filc;Bvol ,Benp ,Ptrfin ,Etp ,Enf ( * ) ,Fr. ,T 1 ,T 2 ,Pv 3 p 1 

DDTMT D>,^1 D*^r, p 4- ^ D + r^ D«« T1 TD Dv,»r^1 

M D V T T 



3904 

3000 

3907 

3909 

■7.a{?iQ 

-/ w/ V ^ 

3917 

3027 

3037 



A0CTCM @Fiie TO 
SUBEND 



Cl ID D, 

w/ W Ly I W I ^ 



BEEP 

TMDIIT "CMTCD CTIC MAMC Tn DC nClCTCn" C,l«« 

Ol'M Wl IW l^L- ,1 

PURGE FileS 



W I « w 

CHTO ■ T.QOtD 

W W ( W -/ U W W 

SUBEND 
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